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19.  ABSTRACT  (Continued) 

Based  on  the  findings  of  this  program  it  is  recommended  that  the  optical  probe  be 
characterized  by  direct  calibration  rather  than  calculated  from  particle  size  distribution 
information. 

A  prototype  measurement  station  was  designed  to  withstand  ideal  and  precursed  shocks 
with  peak  overpressures  up  to  100  psi.  This  preliminary  design  allows  measurements  to 
be  made  at  six  heights  above  the  ground  (V ,  2",  3V,  5",  and  8"). 
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SUMMARY 


l 


Dynamics  Technology's  development  program  has  resulted  in  the  construction  and 
testing  of  a  laboratory  prototype  optical  probe  which  provides  a  measure  of 
the  mass  concentration  of  airborne  dust  in  the  range  of  10~4  to  10“2  gm/cm3 
(p^ust/Pair  ~  t0  !0) -  The  probe  uses  an  infrared  LED  light  source  and  a 
silicon  photodiode  to  measure  the  attenuation  of  light  passing  through  the 
dust.  Optical  fibers  carry  the  light  to  the  measurement  points.  The  use  of 
inexpensive  LEDs,  fibers,  couplers  and  photodiodes  allows  a  low  cost  per  spa¬ 
tial  measurement  point,  and  the  use  of  optical  fibers  allows  delicate  elec¬ 
tronics  and  most  optical  components  to  be  housed  in  a  remote  location.  The 
probe  provides  a  pulsed  optical  signal  allowing  spurious  ambient  light  and 
electronics  zero-drift  to  be  measured  and  accounted  for. 

The  prototype  probe  was  calibrated  over  a  range  of  dust  densities  similar  to 
those  measured  in  the  Pre-Direct  Course  event,  using  both  commercially  avail¬ 
able  graded  sand  and  a  representative  dust  sample  from  the  Direct  Course  test 
bed.  Experiments  demonstrated  that  the  range  of  dust  densities  measured  can 
be  varied  by  changing  the  optical  path  length.  In  addition,  the  particle  den¬ 
sity  and  size  distribution  of  the  graded  sand  and  the  dust  sample  from  the 
Pre-Direct  Course  test  bed  were  measured  by  an  independent  service  laboratory, 
and  the  results  were  used  as  a  basis  for  making  a  priori  estimates  of  the 
optical  attenuation  versus  dust  density. 

A  measurement  station  to  house  the  optical  probes  in  high  explosive  field 
tests  and  in  shock  tube  tests  was  designed.  The  measurement  station  was  sized 
to  minimize  aerodynamic  interference  and  to  withstand  expected  shock  condi¬ 
tions.  A  digital  data  acquisition  system  for  the  probes  was  designed  for  use 
during  high  explosive  field  tests,  as  was  an  underground  vault  for  containing 
the  probe  electronics,  the  optical  emitters  and  detectors,  and  the  data  acqui¬ 
sition  system. 
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This  program  was  conducted  for  the  Defense  Nuclear  Agency  under  contract  DNA- 
000-83-C-0403  to  develop  the  design  for  a  fiber  optic  dust  probe.  The  effort 
was  funded  under  the  Small  Business  Innovative  Research  Program  and  was 
performed  from  July  1983  to  February  1984.  Dr.  Paul  R.  Rohr  was  the  technical 
moni  tor. 

The  Dynamics  Technology,  Inc.  Program  Manager  was  Duane  Hove.  Michael  Petach 
assembled  the  test  equipment,  conducted  the  development  tests  and  prepared  the 
final  designs.  Dr.  C.  Michael  Dube  provided  valuable  technical  guidance. 


This  work  has  undergone  an  extensive  internal  review  before  release,  both  for 
technical  and  non-technical  content,  by  an  independent  review  committee. 
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SECTION  1.  INTRODUCTION 


Dynamics  Technology  has  conducted  a  program  to  develop  a  low  cost  instrument 
for  measuring  airborne  dust  densities.  The  dust  diagnostic  technique  measures 
the  attenuation  of  a  light  beam  across  a  short  pathlength  through  the  dust.  A 
light  source  coupled  to  an  optical  fiber  transmits  the  beam  to  the  measurement 
point,  then  another  fiber  returns  the  attenuated  signal  to  a  photodetector. 

1.1  Objectives 

The  overall  objective  of  this  program  was  to  develop  a  rugged  inexpensive 
optical  probe  to  provide  spatial  (multiple  probe  heights  and  probe  locations) 
and  temporal  measurements  of  airborne  dust  lofted  during  nuclear  blast  wave 
simulations.  Specifically,  there  were  two  major  technical  objectives.  First, 
suitable  components  were  to  be  selected  through  a  laboratory  evaluation  pro¬ 
gram,  assembled  as  a  laboratory  version  of  the  probe,  and  tested  to  determine 
measurement  range  and  accuracy.  Second,  a  preliminary  design  of  a  prototype 
measurement  station  to  house  the  probes  was  to  be  prepared,  and  a  data  acqui¬ 
sition  system  for  field  tests  was  to  be  specified. 


1.2  Background 


During  tests  in  the  1950's  it  was  observed  that  large  quantities  of  dust  were 
lofted  by  the  blast  waves  from  nuclear  bursts.  Extensive  efforts  were  made  to 
quantify  the  effects  of  dust  on  the  blast  flowfield  because  of  the  observed 
enhanced  damage  to  target  structures.  Interest  in  blast  wave  lofted  dust  con¬ 
tinued  because  it  was  found  to  (1)  enhance  the  airblast  induced  dynamic  pres¬ 
sure;  (2)  influence  the  design  of  air  entrainment  systems  for  hardened  shel¬ 
ters;  (3)  obscure  ground  based  radar  and  optical  sensors;  and  (4)  act  as  the 
source  of  dust  for  nonsurface  burst  dust  clouds,  which  can  interfere  with  air¬ 
borne  operations. 


In  order  to  quantify  these  effects,  physical  measurements  of  the  lofted  dust 
from  high  explosive  tests  are  required  to  provide  inputs  to  the  dust  lofting 
models  and  to  check  the  ability  of  predictive  codes  to  calculate  the  dust 
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transport.  Such  measurements  are  difficult  and  the  instruments  are  often  ex¬ 
pensive  because  of  the  high  overpressures  of  interest  (tens  of  atmospheres), 
the  short  durations  of  the  events  (tens  to  hundreds  of  milliseconds),  and  the 
high  dust  masses  involved  (order  of  the  air  density). 

Many  techniques  have  been  considered  to  measure  the  dust  raised  by  the  shock 
wave  during  a  nuclear  blast  simulation.  Electromagnetic  propagation,  ionizing 
radiation,  sound  speed,  pitot  probe  anl  capacitance  gage  techniques  were  cri¬ 
tically  reviewed  by  Hartenbaum.1  Additional  suggestions  for  optical  extinc¬ 
tion,  holography,  microwave  and  x-ray  measurements  have  been  made  by  various 
investigators. 

The  Defense  Nuclear  Agency  conducted  non-nuclear  burst  tests  in  1971  (Middle 
Gust  II)  and  in  1982  and  1983  (Direct  Course  series)  employing  several  dust 
diagnostic  techniques.  The  various  dust  diagnostic  techniques  displayed  vary¬ 
ing  degrees  of  success,  but  none  provided  a  definitive  spatial  and  temporal 
history  of  the  dust  flow.  Two  optical  extinction  dust  measurement  systems 
were  fielded  in  the  Pre  Direct  Course  event;  a  laser  diode-based  optical 
absorption  instrument,2  and  a  HeNe  laser  based  optical  absorption  and  scatter¬ 
ing  system.3  The  laser  diode  system  survived  the  blast  and  provided  a  contin¬ 
uous  history  of  the  dust  flow  past  a  single  point;  the  HeNe  laser  system  only 
provided  an  incomplete  record  because  higher  than  expected  shock  overpressures 
caused  optical  misalignment.  Thus,  the  viability  of  optical  extinction  dust 
diagnostic  techniques  was  demonstrated,  but  the  need  for  a  rugged  system  was 
also  clearly  indicated.  Although  the  laser  diode-system  was  successful,  it 
had  three  drawbacks.  First,  there  was  only  one  spatial  measurement  point. 
Second,  the  measurement  point  was  between  two  blast  wings  such  that  there  is  a 
good  possibility  that  the  measured  flow  was  disturbed.  Third,  this  data  was 
reduced  assuming  a  linear  relationship  between  light  extinction  and  dust  den¬ 
sity  rather  than  the  actual  exponential  relationship.  So  far,  no  single 
approach  has  proven  completely  satisfactory. 


SECTION  2.  THEORY 


Mass  concentration  of  airborne  dust  can  be  determined  by  measuring  the  atten¬ 
uation  of  a  light  beam  that  passes  through  the  dust.  Desired  ranges  of  mass 
concentrations  are  selected  by  varying  the  optical  path  length  through  the 
dust. 


Optical  Transmission  Through  Airborne  Dust 


The  transmission  of  light  through  a  homogeneous  dispersion  of  uniform  diameter 
spherical  particles  which  both  scatter  and  absorb  light  is  given  by  the 
Bouguer-Beer-Lambert  law: 


T  =  =  exp(-  ^  S-j  U)  (1) 

0  K 

where  T  =  transmission  ratio 

I  =  intensity  of  received  light 

I0  =  intensity  of  received  light  when  no  particles  are  present 

Q  =  dimensionless  extinction  factor  which  accounts  for 
both  absorption  and  scattering  losses 

p  =  particle  density 

d  =  particle  diameter 

C  =  mass  concentration  of  particles 

i  =  optical  path  length  through  the  particles 


Real  lofted  dust  flows  have  a  distribution  of  particle  sizes  rather  than  a 
uniform  particle  size.  The  Bouguer-Beer-Lambert  law  can  be  applied  to  a 
polydispersion  of  spherical  particles  in  the  form  given  by  Dobbins4 

T  =  =  exp(-  |  U)  (2) 
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where  Q  =  average  extinction  factor  =  E(NQd2A)/£(Nd2A)  (3) 

d 3 2  =  vol ume-to-surface  mean  particle  diameter  =  ^ 

E(Nd3A)/E(Nd2A)  =  E(xA)/(xA/d) 

and  A  =  particle  size  interval 

d  =  midpoint  of  the  interval  A 
N  =  number  of  particles  per  interval  A 
X  =  %  of  the  total  weight  in  the  interval  A 

The  mass  concentration  as  a  function  of  received  intensity  is  found  by  solving 
equation  (2)  for  C,  yielding 

c  ’  (- 1 pfe  *rl  *"  (5) 

Then  for  a  particle  distribution  characterized  by  the  term  -5 —  the  range  of 

Pd32 

mass  concentrations  determined  by  a  given  range  of  optical  transmissions  can 
be  varied  by  changing  the  interrogating  path  length  (Figure  2.1). 


Once  the  term  — ^ —  is  known,  the  transmission  of  light  across  an  optical  path 
Pd32  3  0 

length  Z  provides  a  measure  of  dust  mass  concentration.  The  term  — ^ —  is  a 

Pd32 

function  of  both  the  particle  size  distribution  and  of  the  optical  character¬ 
istics  of  the  interrogating  instrument.  Obviously,  the  particle  size  distri¬ 
bution  affects  both  Q  and  d32,  as  shown  by  equations  (3)  and  (4).  In  addi¬ 
tion,  the  optical  characteristics  of  the  interrogating  instrument  also  affect 
the  term  Q  in  equation  (3).  For  particle  diameters  larger  than  the  optical 
wavelength  \  (in  the  application  of  interest,  d  »  \),  the  value  of  Q  asymp¬ 
totically  approaches  a  value  of  2  if  no  scattered  light  strikes  the  detector. 
But  0=2  only  if  the  interrogating  beam  is  perfectly  collimated  and  an  ideal 
far  field  detector  is  used.  In  real  devices  the  finite  aperture  of  the  detec¬ 
tor  allows  some  of  the  scattered  light  to  be  measured,  particularly  from  large 
particles  which  scatter  strongly  in  the  forward  direction,  so  Q  may  be  less 
than  2.  Furthermore,  real  dust  particles  are  not  spherical,  so  equation  (2) 
does  not  hold  exactly. 

There  are  several  methods  of  determining  an  optical  probe's  characteristic 
transmission  vs.  dust  mass  concentration  curve  for  a  given  particle  distribu¬ 
tion.  The  most  direct  method  is  to  experimentally  determine  the  term  — 
through  a  calibration  in  which  I/IQ  vs.  C  is  measured  using  a  representative 
dust  sample.  Another  method  is  to  experimentally  determine  d32  and  p  from  a 
representative  dust  sample,  and  assume  a  value  of  0  (typically  Q  »  2  for  large 
particles).  Both  of  these  methods  were  assessed  in  the  present  program. 

An  alternative  approach  to  probe  characterization  is  to  assume  a  value  of  Q> 
measure  p  from  a  representative  soil  sample,  and  then  estimate  d32  from  real 
time  optical  scattering  measurements  made  during  the  event.  This  method  does 
not  depend  on  a  representative  sample  distribution  so  it  has  the  potential  for 
being  more  accurate.  However,  it  is  very  difficult  to  implement,  as  will  be 
discussed  in  the  next  section. 


2.2  In-Si tu  Particle  Sizing  Considerations 


ISM 


In-situ  particle  sizing  can  be  accomplished  by  measuring  the  angular  depen¬ 
dence  of  the  relative  intensity  of  forward  scattered  light.  This  is  possible 
because  forward  scattering  intensity  from  a  particle  of  diameter  d  peaks  at  a 
characteristic  angle  of  0.438  rad.  For  example,  if  near  infrared  light  of  \ 
=  820  nm  is  used,  the  characteristic  scattering  angles  for  particles  in  the 
size  range  of  interest  are  given  in  Table  2.1. 

Table  2.1  Principal  Scattering  Angles  for  Dust  Particles 


PARTICLE  DIAMETER 
d  (microns) 


SCATTERING  ANGLE 
6  (degrees) 


If  a  photodetector  array  located  an  appropriate  distance  (the  order  of  10  cm) 
from  the  scattering  volume  is  used  to  size  10  p  to  100  p  particles  into  4 
ranges,  then  detectors  would  need  to  be  placed  at  0.36  cm,  .14  cm,  .072  cm, 
.048  cm  and  .036  cm  from  the  beam  axis.  Obviously,  the  physical  size  of  com¬ 
mercially  available  photodetectors  makes  this  impractical.  In  addition,  the 
scattered  light  power  levels  are  approximately  20  dB  down  from  the  direct 
beam.  Wuerker2  has  investigated  an  approach  to  scattering  measurements  which 
could  be  adapted  to  the  probe  configuration  developed  in  this  program. 
Wuerker's  system  is  based  on  the  fact  that  all  light  entering  a  lens  at  a 
given  angle  is  brought  to  a  common  focal  point  in  the  focal  plane.  Thus,  the 
conical  scattering  at  an  angle  <t»  is  transformed  by  a  lens  of  focal  length  f 
into  a  ring  of  radius  s  =  f  tan(<t>)  in  the  focal  plane  of  the  lens.  If  the 


focal  plane  is  divided  into  annular  regions,  the  intensity  in  each  zone  is 
proportional  to  the  number  of  particles  in  the  size  range  with  the  correspond¬ 
ing  conical  scattering  angle,  and  also  to  the  cross  sectional  area  of  the  par¬ 
ticles.  However,  the  radii  s  of  interest  for  practical  systems  are  still  the 
order  of  a  millimeter  or  less.  Calculations  in  Appendix  C  indicate  the 
structural  vibration  from  shock  loading  could  result  in  deflections  that  pre¬ 
clude  useful  size  distribution  estimates. 


The  relative  error  in  calculating  the  mass  concentration  C  from  optical  trans¬ 
missions  is 


Thus,  the  relative  error  in  C  is  directly  proportional  to  the  relative  error 
i'n  P,  d32»  i  and  Q.  However,  the  relative  error  is  a  nonlinear  function  of 
intensity.  This  term  can  be  examined  by  plotting  lines  of  constant  measure¬ 
ment  accuracy  versus  and  (Figure  2.2). 

^  *o  o 

This  figure  shows,  for  example,  that  if  a  sensor  has  1%  transmission  resolu¬ 
tion,  then  dust  concentration  measurements  with  10%  accuracy  are  limited  to 
transmissions  ranging  from  90%  to  3%.  Or,  if  a  sensor  measures  transmissions 
down  to  I/I0  =  1%,  dust  concentration  measurements  with  10%  accuracy  require 
at  least  a  1/2%  transmission  resolution. 


SECTION  3.  PROBE  DEVELOPMENT  AND  EVALUATION  PROGRAM 


In  addition  to  cost  considerations,  two  major  design  issues  were  addressed. 
First,  the  emitter  optical  power,  system  power  losses  and  detector  sensitivity 
must  provide  sufficient  dust  density  dynamic  range.  Second,  the  system  must 
physically  fit  into  the  proposed  measurement  station.  The  power  issue  affec¬ 
ted  every  component  in  the  system;  the  problem  of  fitting  the  probe  into  the 
measurement  station  mainly  affected  the  collimation  and  focusing  optics. 


3.1  Design  Considerations 

The  primary  design  criterion  was  that  the  probe  had  to  measure  the  expected 
range  of  dust  mass  concentrations  with  a  reasonable  degree  of  accuracy.  Data 
from  the  Pre-Direct  Course  event  provided  an  estimate  for  the  expected  range 
of  dust  densities. 

For  example,  Wisotski's3  laser-diode  attenuation  data  shows  (once  the  dust 
density  is  calculated  on  an  exponential  scale)  that  normalized  dust  mass  con¬ 
centrations  from  Pdust/Pair  =  t0  pdust/pair  =  9  were  observed  at  a  S1x 
inch  height  above  a  finished  dusty  surface  205  feet  from  ground  zero  (Figure 
3.1).  Dust  densities  decrease  with  distance  above  the  ground,  and  depend  on 
the  distance  from  ground  zero,  so  a  nominal  design  range  of  C  =  10~4  to  10”2 
g/cc  was  chosen.  The  optical  requirements  then  can  be  estimated  following  the 
discussion  in  the  sensitivity  analysis  (Section  2.3).  For  example,  with  an 
instrument  capable  of  1%  optical  resolution,  the  attenuation  caused  by  C  = 
10~4  g/cc  has  to  cause  an  attenuation  of  at  least  5%  in  order  to  achieve  20% 
accuracy.  The  transmission  equation  then  becomes 

-j—  =  0.95  =  exp(a(10-4  g/cc))  (7) 

o 

for  which  a  =  -513  cc/g  (which  can  be  achieved  by  tailoring  the  pathlength, 


TEST  PRf  DIRFCT  COURSE 
EXP  OR I 

STATION  222-  50-LD 


TEST  PRE-DIRECT  COURSE 
EXP  DRI 

STATION  222-.58-L0 


Figure  3.1  Laser-diode  Dust  Measurements 

from  the  Pre-Direct  Course  Event 
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given  the  term  — ^ — ).  A  mass  density  C  =  10-2  g/cc  would  then  cause  a  trans- 
mission  ratio  I / 1 Q  =  6  x  10  3,  which  corresponds  to  a  22  dB  optical  dynamic 
range.  The  necessary  frequency  response  can  also  be  estimated  from  the  Pre- 
Direct  Course  data  of  Figure  3.1,  where  it  is  apparent  that  no  sigificant 
fluctuation  in  dust  density  occurred  in  less  than  0.5  ms.  Thus,  a  2  kHz  samp¬ 
ling  rate  was  considered  adequate.  An  additional  design  goal  was  to  have  the 
capability  of  measuring  and  accounting  for  the  effects  of  non-incident  light 
(e.g.,  scattered  sunlight,  the  fire  ball,  or  glowing  debris)  and  electronic 
zero  drift.  This  can  be  achieved  by  operating  the  probe  with  a  pulsed  optical 
signal . 

The  physical  requirement  for  an  aerodynamic  measurement  station  design  sugges¬ 
ted  using  a  single,  small  diameter  lens  to  collimate  and  to  focus  the  beam. 


3.2  Initial  Component  Selection 


A  literature  survey  showed  that  optical  extinction  measurements  of  dust  densi¬ 
ty  had  been  successfully  achieved  using  a  HeNe  laser,5  a  laser  diode,2  and  an 
infrared  LED6’7  light  sources,  and  photodi odes2 *5 ,6 ,7  and  phototransistors7  as 
detectors.  There  was  no  inherent  performance  advantage  to  any  of  these  sys¬ 
tems  and  none  of  the  systems  incorporated  fiber  optics,  so  it  was  not  evident 
which  type  of  source  or  detector  was  most  suitable.  Therefore,  at  least  one 
of  each  type  of  light  source  and  detector  was  tested  in  the  development  of  the 
probe.  A  variety  of  optical  fibers  with  different  core  diameters  and  core 
materials  were  tested,  since  it  was  uncertain  how  these  parameters  would 
affect  power  transmission  or  coupling  efficiency.  AMP  Optimate  fiber  optic 
connectors  and  device  mounts  were  chosen  since  they  accommodated  the  broadest 
variety  of  fibers,  and  they  were  readily  available  at  low  cost.  The  remaining 
optical  support  equipment  such  as  lenses,  holders  and  positioners  was  avail¬ 
able  in-house.  The  initial  selection  of  components  is  shown  in  Table  3.1. 


DEVICE 


DESCRIPTION 


Sources 

"  '  ‘  ■  * 

— 

-  - - 

1. 

HeUe  Laser 

0.5  raW  X 

*  633  nm.  1  mm  Collimated  Beam 

2. 

Honeywell  SPX  4689-002  1 

Honeywell  SPX  4689-003  t 

In  f rared 

LED  s ,  \  -  820  nm 

,  Fiber  Optic  Mount 

3. 

GE  F5E1 

Infrared 

LED’ s  X  =  850  nm. 

4. 

NEC  NDL30000 

Laser  Diode,  X  -  780  nm  P 

.  =  8.0  mW 

5. 

Mitsubishi  I1L3131 

Laser  Diode,  X  =  820  nm 

Detectors 

1. 

UDT  125  DP/L 

Si  1  icon 

Photodiode  Lensed 

PV 

2. 

UDT  040A 

Sil  icon 

Photodiode  Lensed 

PV 

3. 

GE  F5D1 

Photodarl ington  Fiber  Optic  Case  PV  1 

4. 

Honeywell  SPX  4691-002 

Sil  i con 

Photodiode  Fiber 

Optic  Case  PC 

5. 

Honeywell  SPX  4692-002 

Sil icon 

Photodiode  Fiber  Optic  Case  PC 

6. 

Honeywell  SPX  4690-002 

Integrated  Detector 

i 

Optical  Fibers 

Core  Dia 

.  Attenuation  0 

HA 

I. 

Newport  Research  FML0 

100  u 

4.5  0  850  nm 

.3 

2. 

DuPont  Crofon  0E1040 

lOOOu 

500  0  620  nm 

.53 

3. 

Nissho  Iwai  SH4001 

lOOOu 

300  0  620  nm 

.50 

4. 

General  Fiber  Optics  PS-400-A 

400u 

8  dB  0  820  nm 

.3 

5. 

Quartz  Products  QSF  600A 

600  u 

5  dB  0  820  nm 

.27 

Fiher  Optic  Connectors  and  Couplers 

1. 

AMF  ■ -530530-1 

400  u 

Ferrule 

2. 

AMP  -.30530-5 

600u 

Ferrule 

3. 

AMP  227205-7 

100  u 

Ferrule 

4. 

AMP  ' 30530-2 

lOOOu 

Ferrule 

Devi  _•  Mounts 

I. 

AMP  5  30563-1 

Mount  for  TO- 18,  TO-46,  TO- 

52  Devices 

2. 

AMP  5 '0-565-1 

Mount  for  T1  3/4  Devices 

3. 

AMP  22  7324-1 

Free  Hanging  Coupling 

3.3  General  Development  Plan 

The  selection  of  an  effective  combination  of  optical  components  was  achieved 
through  a  laboratory  evaluation  of  candidate  configurations.  The  first  step 
in  the  evaluation  of  a  configuration  was  to  check  that  sufficient  optical 
power  was  transmitted  through  the  system  and  that  an  adequate  degree  of  col- 
limation  of  the  interrogating  beam  was  achieved.  Then  the  candidate  configu¬ 
rations  were  tested  to  determine  their  extinction  coefficient  and  their  per¬ 
formance  characteri sti cs.  Specifically,  the  transmission  ratio  I / I Q  across  a 
known  path  length  was  measured  at  several  independently  measured  dust  mass 
concentrations.  This  data  was  plotted  in  the  form  I / 1 0  =  exp  (aC),  and  a 
least  squares  fit  of  the  transmission  equation  (5)  was  performed. 


The  path  length  l  could  be  factored  out  of  this  curve  fit,  allowing  the 
probe's  measurement  range  to  be  predicted  as  a  function  of  path  length.  This 
path  length  scaling  was  verified  by  performing  calibrations  using  different 
path  lengths.  In  addition,  the  dust  sample's  particle  distribution  and  speci¬ 
fic  gravity  were  measured  by  an  independent  service  laboratory,  and  the  terms 
i,  p  and  d32  were  used  to  provide  an  a  priori  estimate  of  a  for  comparison 
with  the  measured  laboratory  values. 


3.4  Test  Facility  and  Calibration  Procedure 

A  test  facility  which  provided  a  controlled  density  dust  flow  was  designed  and 
fabricated.  This  facility's  design  was  based  upon  the  dust  flow  calibration 
apparatus  used  by  Ausherman.5  In  brief,  dust  flows  down  a  vertical  channel, 
past  the  probe  and  through  a  tube  of  known  volume.  The  dust  density  is  mea¬ 
sured  by  simultaneously  closing  off  the  top  and  bottom  of  this  tube,  weighing 
the  captured  dust  and  dividing  the  dust  weight  by  the  tube  volume. 

The  test  facility  (Figure  3.2)  consists  of  an  eight  foot  high  vertical  flow 
channel  with  an  8  inch  by  8  inch  square  cross  section.  The  front  face  of  the 

flow  channel  is  plexiglass,  allowing  the  dust  flow  to  be  observed.  A  1,000  in3 

capacity  (~  50  lbs)  hopper  at  the  top  of  the  flow  channel  holds  the  dust 

sample.  The  dust  flow  out  of  the  hopper  is  controlled  by  a  variable  orifice 
made  of  two  perforated  plates  which  slide  relative  to  each  other.  As  the 
plates  slide,  the  holes  on  the  two  plates  overlap,  allowing  dust  to  flow  out 
in  a  uniform  distribution  across  the  entire  flow  channel  cross  section. 

Adjustable  stops  on  the  plates  control  the  degree  of  overlap,  and  thus  the 

flow  rate.  Dust  mass  concentrations  from  ~  2  x  10-4  to  8  x  10~3  g/cc  were 
achieved  in  this  manner.  As  the  dust  falls  down  the  flow  channel  it  is  uni¬ 
formly  dispersed  by  screens  at  4  inch  intervals.  These  screens  are  removable 
so  that  different  mesh  sizes  can  be  used  with  different  soil  types.  The 
optical  probe  protrudes  into  the  flow  through  slots  in  the  flow  channel  wall 
located  just  below  the  lowest  screen.  The  optical  probe  is  aligned  and  held 
at  a  preset  spacing  by  magnetic  bases  on  a  metal  table  surface  outside  the 
flow  channel  (Figure  3.3). 
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Figure  3.2  Dust  Flow  Test  Facility 


The  dust  is  collected  in  an  acrylic  tube  (5.75  inches  in  diameter  by  23.125 
inches  tall)  held  in  a  vertical  orientation  concentric  with  the  flow  channel 
walls.  Dust  flows  freely  around  and  through  the  tube  until  a  dust  density 
measurement  is  made  by  simultaneously  closing  off  the  top  and  bottom  of  the 
tube  with  sliding  blades.  The  dust  flowing  through  the  tube  prior  to  closure 
is  trapped  and  falls  into  a  container  held  on  the  bottom  blade. 

The  test  procedure  used  to  acquire  light  transmission  vs.  dust  density  data 
was  as  follows.  First  the  stops  on  the  hopper  orifice  were  set  for  an  appro¬ 
priate  flow  rate.  Before  initiating  the  dust  flow,  the  detector  electronics 
were  zeroed  (if  necessary),  the  optical  windows  were  cleaned,  and  the  voltage 
VQ  at  the  unattenuated  intensity  IQ  was  recorded. 


*>« 
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Figure  3.3  Dust  Probe  Placement  in  Test  Facility 

After  the  hopper  orifice  opened,  dust  was  allowed  to  flow  freely  until  the 
flow  appeared  well  dispersed,  and  the  probe  voltage  output  displayed  on  a 
digital  multimeter  had  stabilized.  This  usually  only  took  a  few  seconds.  The 
computer  controlled  data  acquisition  program  then  digitally  sampled  the  probe 
output  voltage  (typically  600  samples  at  300  Hz).  When  the  data  acquisition 
was  completed,  the  dust  collection  tube  was  immediately  sealed  off.  The  data 
acquisition  time  span  was  set  so  that  a  particle  falling  at  slightly  less  than 
the  terminal  velocity  of  a  sphere  with  the  mean  dust  particle  diameter  (30 
cm/sec)  would  travel  the  height  of  the  dust  collector  tube.  This  meant  that  a 
particle  passing  the  probe  at  the  beginning  of  the  data  acquisition  would  just 
reach  the  bottom  of  the  tube  at  the  end  of  the  data  acquisition,  and  the  dust 
trapped  in  the  tube  would  correspond  with  the  dust  that  had  contributed  to  the 
transmission  measurement. 


The  post-flow  value  of  VQ  was  also  recorded  to  assess  dust  adhesion  to  the 
windows.  If  there  was  a  significant  decrease  in  VQ  due  to  dust  coating,  the 
post-flow  value  of  VQ  was  used.  The  value  of  VQ  was  entered  into  the  compu¬ 
ter,  which  then  plotted  the  normalized  voltage  history  V/VQ  =  I/I  and  compu¬ 
ted  the  average  and  rms  normalized  voltages.  Typical  hard  copy  outputs  are 
shown  in  Figure  3.4. 


16 


6 
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TIME  (sec  > 

DIODE  OUTPUT  US  TIME 


PUN  NUMBER  2? 

MEAN  I'lo  -  +2. 2441E-002 

Cl  =  ♦ 1 . 0000E+000 

PMS  I/Io  «  ♦2.5128E-003 

C2  *  8 . 0000E+000 

MEASURED  WEIGHT  =  70.01 

FREQUENCY  =*  300 

CALCULATED  DENSITY  •  +7. 1970E-003 

•  0F  SAMPLES  -  601 

MEAN  UOLTAGE  -  +2.2441E-00I 

TPUT  US  TIME 


RUN  NUMBER  6 

MEAN  UIo  »  +3.  1631E-001 

Cl  «=  ♦1.0006E+000 

RMS  Ulo  ‘  *  1 . 5743E-002 

C2  »  0 . 0000E *300 

MEASURED  WEIGHT  -  10.24 

FREQUENCY  «  300 

CALCULATED  DENSITY  ■  +1.0327E-003 

•  0F  SAMPLES  -  601 

MEAN  UOLTAGE  ■  ♦3.1681E*000 

fi] 


SECTION  4.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


During  the  experimental  development  of  the  optical  dust  probe,  several  com¬ 
binations  of  optical  emitters,  optical  fibers,  collimating  optics,  and  detec¬ 
tors  were  tested.  These  tests  and  the  pertinent  results  are  summarized  in 
Table  4.1 

The  dust  density  vs.  transmission  test  results  obtained  with  the  preliminary 
probe  configurations  are  shown  in  Figures  4.1  through  4.3  (the  final  configu¬ 
ration  will  be  discussed  later).  It  is  apparent  that  each  configuration  pro¬ 
vided  data  that  followed  the  form  of  the  transmission  equation  I / 1 Q  =  e“a^. 

Figure  4.3  shows  that  the  calibration  curves  of  a  given  probe  configuration 
scale  with  a  change  in  the  optical  path  length,  as  expected. 

During  the  development  program,  the  evaluation  of  the  various  components  pro¬ 
vided  some  clear  guidelines  for  the  selection  of  the  final  probe  configura¬ 
tion.  For  example,  it  was  found  that  while  both  a  HeNe  laser  and  an  infrared 
LED  provided  sufficient  optical  dynamic  range  to  measure  the  desired  range  of 
dust  densities,  the  LED  was  clearly  superior  for  the  present  application  since 
i  t 


(1)  is  less  expensive  (~  $20  vs.  ~  $370), 

(2)  is  much  more  compact, 

(3)  requires  less  power,  e.g.,  can  be  powered  by  batteries 

(4)  is  more  easily  coupled  to  an  optical  fiber  (the  chosen  LED  is 
factory  mounted  in  a  fiber  compatible  housing) 

(5)  can  be  electronically  pulsed  to  allow  ambient  light  noise  and 
electronic  zero  drift  to  be  measured. 

A  GaAlAs  infrared  (820  qm)  LED  was  chosen  as  the  optical  emitter  (Figure  4.4). 
It  was  driven  by  a  simple  transistor  switching  circuit  which  allowed  contin¬ 
uous  or  pulsed  operation  (Figure  4.5). 
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IR  LED  Ferruled  Collimating  lens  Ferruled  UDT  125  DP/L  Demonstrated  pulsed  operation 

Honeywell  General  Fiber  Optics  17.6  mm  dla.  15.9  mm  f.l.  Quartz  FVoducts  Silicon  Riotodlode  Illustrated  optical  dynamic 

SPX  4689-003  PS-400-A  Focusing  lens  QSF-600A  In  Trans  Impedance  range  and  system  frequency 

Msed  Operation  18  mm  dla.  33  mm  f.l.  Circuit  response. 
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4.3  Experimental  Demonstration  of  Path  Length  Scaling 
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Figure  4.5  LED  Driver  Circuitry 


Experiments  with  coupling  different  fibers  to  the  LED  showed  that  there  was  a 
trade-off  between  the  power  launched  through  an  optical  fiber  of  a  given  core 
diameter  and  the  degree  of  collimation  attainable  with  the  fiber.  This  is 
because  the  larger  the  fiber  diameter,  the  less  it  acts  as  a  point  source.  A 
400  pm  diameter  transmitting  fiber  provided  a  good  compromise.  A  600  core 
was  used  for  the  receiving  fiber  since  the  larger  fiber  face  reduced  vibration 
sensitivity  and  increased  coupling  efficiency.  Both  optical  fibers  were  term¬ 
inated  with  AMP  Optimate  ferrules,  allowing  them  to  be  easily  and  repea tably 
coupled  (Figure  4.6) . 
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Figure  4.6  Fiber  Termination  and  Connector 


A  silicon  photodiode  operating  in  the  photovoltaic  mode  in  a  transimpedance 
circuit  was  chosen  as  the  optical  detector.  The  output  voltage  of  this  circuit 
is  linearly  proportional  to  intensity  (Appendix  C)  which  allows  direct  deter¬ 
mination  of  the  transmission  ratio  from  measured  output  voltage  ratio.  The 
p  todiode  was  mounted  in  an  AMP  active  device  mount  (Figure  4.7)  and  wired 
directly  into  a  transimpedance  amplifier  circuit  (Figure  4.3)  to  reduce  elec¬ 
tronic  noise  and  increase  the  frequency  response. 
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The  final  fiber  optic  dust  probe  configuration  was  calibrated  using  graded 
Ottawa  Sand  and  a  soil  sample  from  the  Pre-Direct  Course  test  bed  (Figure 
4.10).  In  each  case,  the  measured  transmission  ratio  v?  measured  dust  mass 
concentration  data  was  fit  (in  a  least  squares  sense)  th  a  curves  of  the 
form 


C  =  An  ( I/I  )  /otJt  .  (a) 

For  the  calibrations  with  Ottawa  Sand,  the  best  fit  is 

C  =  An  (I/I  J/-209.  5  A  .  (9) 

The  relative  deviations  about  the  best  fit  curve  curve  were  calculated,  and  it 
was  found  that  902  of  the  points  fell  within  ±  302  of  the  curve  fit. 

Similarly,  the  best  fit  of  the  calibrations  using  the  Pre-Direct  Course  soil 
sample  are 


C  =  An  { I/I  )/-93. 5  A  (10) 

In  this  ;ase,  902  of  the  data  fell  within  ±252  of  the  cr  ve. 

Thus,  based  on  these  calibrations,  the  natural  logarithms  of  the  transmission 
ratio  as  a  predictor  of  dust  mass  concentration  provides  a  accuracy  of  approx¬ 
imately  ±302. 

Extinction  curves  for  the  Pre  Direct  Course  soil  sample  and  the  Ottawa  sand 
were  also  calculated  from  the  average  of  the  pre-test  and  post-test  particle 
size  distribution  and  particle  specific  gravity  data  returned  by  the  indepen¬ 
dent  testing  laboratory  (Appendix  A).  Putting  these  results  in  the  form  of 
the  extinction  equation  and  assuming  Q  =  2,  the  transmission  equation  for 
Ottawa  sand  yields 

C  =  in  (I/I0)/-H9.9  A  (11) 


i 
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and  the  transmission  equation  for  the  Pre-Direct  Course  soil  sample  yields 

C  =  An  ;i/IQ)/-57.6  A  (12) 

The  transmission  data  for  the  two  dust  types  varies  with  the  term — as 

d3  2P 

expected.  This  is  evident  in  the  fact  that  the  ratio  of  the  extinction  coef¬ 
ficient  terms  for  the  two  dust  samples  as  determined  by  the  soil  analysis  of 

the  term  —  i  s 
032P 

Ottawa  _  119.4  _  9  ,  mm 

Direct  Course  "  57.6  "  ' 

and  the  calibrations  yield  the  extinction  coefficient  ratio 


Ottawa  _  208.7 

Direct  Course  92 . 6 


However,  it  is  evident  that  the  extinction  coefficients  calculated  from  soil 
analysis  differed  markedly  from  the  extinction  coefficients  determined  by 
calibration,  as  shown  in  Figure  4.11.  For  the  Ottawa  sand,  the  dust  mass 
concentration  at  a  given  transmission  as  calculated  from  the  soil  analysis  is 
75*  higher  than  the  dust  mass  concentration  determined  by  the  calibration 
curve.  For  the  Direct  Course  dust  sample,  the  dust  mass  concentration  calcu¬ 
lated  from  the  soil  analysis  is  62%  higher  than  the  dust  mass  concentration 
determined  by  the  calibration  curve. 

These  differences  are  not  explained  by  the  variation  of  0  between  1  and  2 
since  varying  the  value  of  Q  from  the  assumed  value  of  2  to  1  moves  the  calcu¬ 
lated  values  of  the  extinction  coefficient  further  away  from  the  measured 
values.  Nor  can  they  be  satisfactorily  explained  by  the  uncertainties  in  the 
two  different  techniques,  as  discussed  in  Appendix  B.  The  exact  cause  of 
these  differences  is  not  presently  known,  however  it  is  interesting  to  note 
that  Cashdollar  et  al.7  found  that  their  dust  probe  calibrations  of  transmis¬ 
sion  vs.  dust  density  using  a  similar  dust  probe  also  fell  below  the  curve 
predicted  by  particle  size  and  density  information  (Figure  4.12).  Based  on 
the  results  of  the  present  program  it  is  clear  that  to  maximize  the  accuracy 
of  the  present  probe  it  must  be  characterized  by  calibration  with  a  represen¬ 
tative  dust  sample. 
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Figure  4.11  Comparison  of  Extinction  Curves  Determined  by  Calibration 
and  Extinction  Curves  Calculated  From  Soil  Analysis 
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The  pulsed  operation  of  the  probe  is  illustrated  in  Figure  4.13,  which  shows  a 
1  kHz  pulsed  optical  signal  transmitted  through  a  dust  flow  of  approximately 
pdust'/oair  =  3.  The  intensity  ratio  indicated  when  the  LED  is  turned  off 
represents  the  total  system  noise  level,  including  ambient  light  and  electron¬ 
ic  noise.  The  measurement  resolution  is  evident  in  the  bit  flip  noise  when 
the  LED  is  off.  This  figure  shows  the  probes  30  dB  dynamic  range,  and  the 
rise  and  fall  times  of  the  system  (approximately  10  n  sec). 
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Figure  4.13  Typical  Pulsed  Operation  Dust  Probe  Output 
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SECTION  5.  MEASUREMENT  STATION  DESIGN 


A  prototype  measurement  station  to  house  the  optical  dust  probes  during  both 
high  explosive  field  tests  and  shock  tube  tests  was  designed.  There  were  two 
major  requirements  that  guided  this  design.  First,  the  measurement  station 
must  be  able  to  withstand  shock  and  dynamic  pressure  loads  without  allowing  a 
degradation  of  probe  performance.  Second,  the  measurement  station  design 
should  minimize  aerodynamic  disturbance  of  the  flowfield  at  the  measurement 
points.  As  an  aide  to  meeting  these  requirements,  the  design  was  based  on  the 
measurement  stations  successfully  fielded  by  Carpenter,  Hove  and  Batt  in  the 
Middle  Gust  III  and  IV  boundary-layer  velocity  profile  tests.9 

Briefly,  the  measurement  station  consists  of  two  steel  blades  which  protrude 
vertically  into  the  flowfield  and  contain  six  optical  dust  probes  at  various 
heights.  One  blade  contains  the  transmitting  optics  (a  single  lens  collima¬ 
tor)  for  all  six  probes  and  the  other  blade  contains  the  receiving  optics  la 
single  lens  collector)  for  all  six  probes.  The  blades  are  bolted  to  a  steel 
base  plate  at  a  separation  which  provides  the  appropriate  optical  path  length 
for  the  expected  dust  densities.  In  shock  tube  tests,  the  optical  path  length 
can  be  modified  as  the  tests  progress  to  optimize  the  probe's  measurement 
range.  A  hollow  steel  sting  extends  horizontally  forward  from  each  blade  into 
the  dusty  flow  at  each  measurement  height.  The  collimated  light  beam  travels 
through  the  sting  to  a  mirrored  steel  insert  which  reflects  the  beam  through 
sapphire  windows,  across  the  dusty  flow  to  the  other  identical  sting. 

In  shock  tube  applications,  the  two  vertical  blades  and  base  plate  would  be 
mounted  on  the  shock  tube  floor.  In  field  applications,  the  blades  and  base 
plate  would  be  bolted  to  a  buried  concrete  pad.  This  concrete  pad  would  con¬ 
tain  the  optical  emitters  and  detectors,  and  a  data  acquisition  system. 
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5.1  Dust  Probe  Rake 


The  dust  probe  rake  blades  are  15  inch  tall  by  12  inch  long  by  2  inch  thick 
steel  plates  with  a  60°  knife  edge  on  the  front  face  (Figure  5.1).  Six  1  inch 
diameter  holes  are  drilled  in  the  front  face  to  hold  the  lens  fixtures  and 
probe  stings,  as  shown  in  Figure  5.2.  The  base  plate  and  blades  are  presently 
designed  to  provide  optical  path  lengths  of  2.5  cm  and  5  cm  (Figure  5.3). 
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Figure  5.1  Dust  Probe  Rake  Blades 


Figure  5.2  Sting  and  Optical  Fixture  Mounted  in  Rake  Blade 

The  sizing  of  the  probe  stings  was  influenced  by  the  diameter  of  the  collima¬ 
ted  beam,  aerodynamic  interference  considerations,  and  structural  require¬ 
ments.  The  sting's  outer  diameter  determined  the  minimum  probe  separation 
(both  vertically  and  horizontally)  due  to  the  aerodynamic  restriction  that  the 
minimum  probe  separation  be  at  least  2  probe  diameters.  The  sting's  length 
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was,  in  turn,  influenced  by  the  diameter  through  the  aerodynamic  rule-of-thumb 
that  the  probe  support  structure  should  be  about  10  probe  diameters  from  the 
probe  tip.  Since  the  expected  optical  path  length  (horizontal  separation)  is 
~  2  inches,  the  maximum  sting  diameter  is  1  inch.  The  10  mm  beam  diameter 
used  in  the  final  optical  probe  configuration  determined  the  inner  diameter  of 
the  probe  sting  tubes.  Based  on  these  considerations  and  the  expected  blast 
loads  (as  shown  in  Appendix  C)  the  stings  shown  in  Figure  5.2  were  specified. 
A  simple  dynamic  structural  analysis  was  performed  to  check  the  maximum  stress 
and  deflections  of  the  stings  and  blades  under  ideal  and  precursor  shock  con- 
di  tions. 


Figure  5.3  Rake  Base  Plate 


5.2  Electronics  and  Instrumentation  Vault 


During  high  explosive  field  tests  the  LEDs,  LED  drivers,  photodiodes,  trans¬ 
impedance  amplifiers,  and  data  acquisition  system  will  be  shock-mounted  in 
underground  steel  vaults  to  protect  them  from  the  blast  and  ground  shock.  The 
vault  design  is  similar  to  that  fielded  by  Carpenter,  Hove  and  Batt8  in  the 
blast  wave  boundary-layer  measurements  for  Middle  Gust  III.  The  vaults 
(Figure  5.4)  are  36  inch  diameter  by  36  inch  deep  steel  cylinders  with  1/4 
inch  thick  walls  and  bottom.  A  2  inch  wide  by  1  inch  thick  flange  welded  to 
the  top  of  the  cylinder  provides  an  anchor  for  a  40  inch  diameter  by  1-1/2 
inch  thick  steel  lid.  The  vaults  will  be  cast  into  the  same  concrete  block 
that  the  probe  rake  is  bolted  to,  with  the  vault  lid  flush  with  the  block 
surface . 


Rough  estimates  of  the  necessary  space  for  the  electronics  and  data  acquisi¬ 
tion  circuitry  are:  one  8"  x  4"  board  to  hold  the  6  LEDs  with  their  driving 
circuitry,  one  8"  x  4"  board  to  hold  the  6  LEDs  with  their  transimpedance 
amplifiers,  and  one  12"  x  12"  board  to  hold  the  digital  data  acquisition  and 
timing  circuits.  If  these  components  are  mounted  in  a  12"  cube  there  is  at 
least  1  foot  of  free  space  around  the  electronics,  which  should  provide  ample 
room  for  batteries  and  shock  mounts. 

5.3  Digital  Data  Acquisition  System 

The  data  acquisition  system  for  field  tests  was  based  on  an  existing  DynaTech 
portable  digital  data  acquisition  system.  In  the  present  design,  the  data 
acquisition  system  is  a  self  contained,  battery  powered  unit.  The  system 
digitally  samples  the  photodiode  voltage  and  stores  the  data  in  dynamic  RAM 
under  the  control  of  a  mi  coprocessor.  After  the  event,  the  data  is  loaded 
onto  a  portable  digital  tape  recorder  and  later  reduced  on  a  minicomputer. 
The  components  in  the  expected  configuration  are  shown  in  Figures  5.5.  An 
Intel  80C31  microprocessor  uses  software  stored  in  A2716  EPROM  to  control  the 
system.  The  data  acquisition  process  is  started  when  an  external  signal  (from 
the  test  operator  or  a  separate  pressure  transducer)  provides  an  interrupt. 
Then  an  Analogic  MP6112-DR  ADC  samples  the  6  dust  probe  signals  at  5  kHz  with 
8  bit  (1/2  %)  resolution.  The  LEDs  are  also  under  the  microprocessor  control 
so  that  the  chopped  light  signal  is  synchronized  with  the  conversion  timing. 
The  data  is  stored  in  64  K  dynamic  RAM  chips.  Each  chip  holds  approximately  2 
seconds  of  data  from  the  six  probes. 

Since  it  is  desirable  to  have  I0  information  from  both  before  and  after  the 
event  to  determine  if  the  optical  windows  become  coated  with  dust,  time 
requirements  dictate  the  use  of  at  least  two  64  K  RAM  chips. 
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SECTION  6.  CONCLUSIONS  AND  RECOMMENDATIONS 


Dynamics  Technology  has  completed  the  phase  one  development  of  a  fiber  optic 
light  extinction  probe  for  measuring  the  mass  concentration  of  dust  lofted  in 
high  explosive  field  tests  and  shock  tube  experiments.  It  was  found  that  low 
cost  optical  components  combined  with  fiber  optics  provided  the  necessary 
optical  dynamic  range  for  measuring  mass  concentrations  in  the  range  of  10~4 
to  10“2  gm/cm3  (pdUSt/Pair  ~  0,1  t0  10^  Calibrations  using  representative 
dust  samples  demonstrated  that  the  probe  could  make  repeatable  dust  mass 
concentration  measurements  in  this  range.  In  these  calibrations,  90*  of  the 
mass  concentration  measurements  fell  within  30%  of  the  calibration  curves. 
Laboratory  tests  demonstrated  that  the  measured  range  of  dust  densities  could 
be  easily  varied  by  changing  the  optical  path  length.  Information  from  an 
independent  service  laboratory  determined  that  the  probe's  calibration  curves 
scaled  to  the  particle  size  distributions  of  the  representative  dust  samples 
as  expected.  However,  estimates  of  dust  mass  concentration  based  on  the 
particle  size  distribution  information  differed  from  estimates  based  on  the 
calibrations  by  up  to  75%.  Therefore,  based  on  the  findings  of  this  program, 
it  is  recommended  that  the  optical  probe  be  characterized  by  direct  calibra¬ 
tion  rather  than  by  particle  size  distribution  information. 

The  optical  dust  pr^,be  in  its  present  configuration  is  ready  for  use  in  dust 
lofting  experiments.  A  minor  hardware  substitution  could  improve  the  probe's 
performance  at  little  additional  cost.  This  would  simply  consist  of  using  a 
Honeywell  SPX  4689-004  LED  instead  of  a  Honeywell  SPX  4689-003  LED.  (The  SPX 
4639-004  was  not  used  for  the  experimental  prototype  was  that  it  was  not 
available  in  time).  This  LED  emits  about  60%  more  optical  power,  and  though 
the  present  probe  demonstrated  ample  dynamic  range,  an  increase  in  dynamic 
range  is  always  desirable. 
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A  prototype  measurement  station  was  designed  which  allows  dust  mass  concentra¬ 
tion  measurements  to  be  made  in  ideal  and  precursed  shocks  with  overpressures 
up  to  100  psi.  At  each  station,  measurements  are  made  at  six  heights  above 
the  ground  (1/2",  2",  3-1/2",  5",  8"  and  12").  The  measurement  station  design 
allows  the  optical  path  length  to  be  changed  between  tests,  thus  changing  the 
measurement  range.  In  addition,  the  design  includes  an  underground  vault  for 
housing  the  opto-electroni cs  during  high  explosive  field  tests,  and  a  self 
contained  data  acquisition  system  for  such  tests. 

The  combination  of  the  fiber  optic  dust  probe  and  the  measurement  station  are 
suitable  for  the  upcoming  CERF  shock  tube  tests  and  the  Miners  Scale  Field 
tests.  It  is  recommended  that  two  stations  be  fielded,  one  at  a  range  corres¬ 
ponding  to  30  psi  overpressure,  the  other  at  a  60  psi  overpressure,  in  both 
these  tests. 
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APPENDIX  A.  SOIL  SAMPLE  ANALYSIS 


Both  representative  dust  samples  were  analyzed  by  an  independent  service 
laboratory  to  determine  their  particle  density  p  and  their  mean  surface 
weighted  diameter  6-32-  The  Pre-Direct  Course  sample  was  analyzed  before  and 
after  the  dust  probe  tests.  The  Ottawa  Sand  was  only  analyzed  after  the  dust 
probe  tests;  prior  to  the  tests  the  suppliers  data  sheet  was  used.  The  par¬ 
ticle  density  p  was  given  by  the  bulk  specific  gravity,  and  the  mean  surface 
weighted  diameter  was  determined  by 


= 


(A .  1 ) 


The  data  and  calculations  follow.  A  summary  of  the  results  is  given  below 


Table  A.l  Summary  of  Soil  Analysis  Results 


PRE-TEST 

POST  TEST 

Ottawa  Sand 

d32  =  0.0078  cm 

d32  =  0.0112  cm 

S. 

G.  =  Not  given 

S.G.  =  2.67 

Pre-Direct  Course  Soil 

1 

2  =  0.0190  cm 

G.  =  2.57 

d32  =  0.0213  cm 

S.G.  =  2.59 

Table  A-2  Independent  Soils  Laboratory  Results 
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Table  A-3.  Mean  Surface  Weighted  Particle  Diameter  Calculations 


INTERVAL  mm 


Pre-DIrect  Course  Sample,  Pre-test 


8 

~ 

.991 

1.381 

1.6855 

1. 

1.381 

.819 

91 

- 

.589 

.402 

.79 

8. 

3.22 

4.071 

89 

- 

.295 

.294 

0.442 

11. 

3.23 

7.317 

95 

- 

.147 

.148 

.221 

17. 

2.52 

11.38 

47 

- 

.074 

.073 

.1105 

21. 

1.53 

13.87 

74 

- 

.040 

.034 

.057 

29.7 

1.01 

17.72 

4 

- 

.026 

.014 

.033 

1.3 

.0182 

0.552 

26 

- 

.015 

0.11 

.0205 

1.1 

.012 

0.590 

15 

- 

.01 

.005 

.0125 

2.5 

.013 

1.000 

1 

- 

.0074 

.0026 

.0087 

5.9 

.0153 

1.76 

074 

- 

.0036 

.0038 

.0055 

0.4 

.00152 

0.276 

036 

- 

.0015 

.0021 

.00255 

0.6 

.00126 

0.494 

015 

- 

.0 

.0015 

.00075 

0.5 

.00075 

1.000 

co  ir> 


Table  A-3. 

Mean  Surface  Weighted  Particle 
Calculations  (continued) 

Diameter 

Ottawa  Sand,  Pre-Test 

INTERVAL  mm 

A  mm  d  mm  x 

_ X  A 

.420 

-  .297 

.123 

.3585 

0.4 

.003 

.008 

.297 

-  .210 

.087 

.2535 

2.0 

.174 

.686 

.210 

-  .149 

.1795 

8.0 

.488 

2.719 

.149 

-  .105 

.044 

.1270 

35.5 

1.562 

12.299 

.105 

-  .074 

.031 

.0895 

31.0 

.961 

10.737 

.074 

-  .053 

.021 

.0635 

12.4 

.2604 

4.100 

.053 

-  .00 

.053 

.0265 

10.5 

.  5565 

21.000 

~Z 

=  4.005 

E  =  51.55 

S.G.  =  Not  Given 


d32  =  0.0078  cm 


Ottawa  Sand,  Post-Test 


INTERVAL  mm 


94 

.442 

1. 

•  ' 

48 

.221 

11. 

1. 

73 

.1105 

57. 

4. 

34 

.057 

30.4 

1. 

.12  E 


32  "  0.0 


APPENDIX  B.  UNCERTAINTY  ANALYSIS 


Uncertainties  in  Extinction  Curves  Obtained  Through  Direct  Calibration 

The  measured  quantities  in  the  direct  calibration  of  the  optical  probe  are: 

(1)  the  transmission  ratio  I/IQ  measured  as  V/VQ  by  a  photodiode 

(2)  the  mass  concentration  C  as  measured  by  the  dust  collection  tube 

(3)  the  path  length  A  is  measured  with  calipers. 

The  measurement  of  I / I 0  is  made  with  a  photodiode  in  a  transimpedance  circuit 
so  the  output  voltage  is  linearly  related  to  the  received  intensity  and  I/IQ  = 
V/VQ.  It  is  important,  then,  to  determine  whether  the  relationship  between  I 
and  V  is  linear  to  within  a  certain  percentage  over  the  measurement  range  of 
interest.  This  can  be  done  by  examining  the  electrical  equivalent  circuit  for 
a  PV  photodiode  (Figure  B.l).  The  flow  of  photocurrent  develops  a  load  volt¬ 
age  which  will  forward  bias  the  ideal  diode,  which  limits  the  PV  photodiode 
response  linearity. 


rloao 


symbols 

Ip  »  PHOTOCURRENT 

Cj  s  JUNCTION  CAPACITY 

In  *  equivalent  noise  current 

RSM  *  SHUNT  RESISTANCE 

Vs  «  SIGNAL  voltage 

VN  *  NOISE  VOLTAGE 

Rj  «  SERIES  RESISTANCE 


Figure  B.l  Photovoltaic  Photodiode  -  Electrical  Model 


B-l 


For  a  given  degree  of  linearity,  the  maximum  photocurrent  is  given  by 


■P.  -  “  10'3  A  (B.D 

where: 

Ipm  =  Maximum  photocurrent  for  P  linearity,  A 

P  =  Percent  linearity  desired  *  1002 

(e.g.,  for  12  linearity,  P  =  0.01) 

Rs  =  Series  Resistance,  Q 

RL  =  Load  Resistance,  Q 

RS[.|  =  Shunt  Resistance,  Q 

The  use  of  an  operational  amplifier  as  a  "current-to-vol tage"  converter, 
(Figure  B.2)  creates  a  unique  solution  to  the  linearity  limitations  imposed  on 
PV  photodiodes  by  the  terminating  load  impedance.  In  this  transimpedance  con¬ 
figuration,  the  PV  photodiode  views  a  load  impedance,  equal  to  the  feedback 
impedance  divided  by  the  open-loop  amplifier  gain. 


Figure  B.2  Photodiode  Transimpedance  Circuit 


(B.2) 


a  '  A{1  T  MfCf)2]1/2 

A  very  good  approximation  of  photodiode/operational  amplifier  response  linear¬ 
ity  can  be  determined  by  the  substitution  of  (Za)  for  RL  in  Equation  (B.l). 


Since  the  voltage  output  of  the  transimpedance  amplifier  is  Vq  =  I^R^.  , 
equation  (1)  can  be  rewritten  to  give  the  maximum  voltage  output  within  a 
given  degree  of  linearity  as 


V 

o 

max 


I  R,  =  R,  [ 
pm  f  f  L 


25  x  10_ 3 1  „„  r  P  X  RSH  i 

J  *n  L—  z(ayJ 


R$  +  Z(a) 


■Rs  + 


(B.3) 


for  typical  component  values. 


A 

=  105 

(B.4) 

RS 

=  1  Q  (lead  resistance) 

(B.  5) 

rsh 

=  40  MQ 

(B.6) 

Rf 

=  500  k Q 

(B.  7) 

cf 

=  2  pf 

(B  .8 ) 

0) 

=  62,831  r/s  (f  =  10  kHz) 

(B.9) 

the  load  impedance  Z(a)  is 
Z(a )  = 


5  x  105  Q 


10b [ 1  +  (6.3  x  104  x  5  x  10b  x  2  x  10_12)2]i/Z 


5.0  Q 


(B. 10) 


and  the  maximum  voltage  output  is 


V  =  5  x  105 
o 

max 


r25  x  10-3 
L  1  +  5  Q  J 


.  rO.Ol  x  5  x  107  Q-i 

[ — rm - ] 


=  36  vol ts 


(B. 11) 


Therefore,  the  typical  photovoltaic  photodiode  in  a  transimpedance  circuit  can 
easily  provide  a  0  -  10V  output  signal  that  provides  a  linearity  within  1% 
over  the  entire  range. 


The  measured  V/VQ  reflects  I / 1 0  to  within  1%  as  verified  with  ND  filters. 
V/VQ  fluctuated  over  the  course  of  the  data  acquisition  period  (Figure  B.l) 
with  the  degree  of  fluctuation  increasing  with  increasing  dust  density.  In 
general,  these  fluctuations  were  on  the  order  of  ±  10%. 


i 

✓ 

I 

0 

1.0E-J 


DIODE  OUTPUT  US  TIKE 


RUN  NUMBER  IS 

MEAN  UIO  -  ♦!.  I376E-BBI 

ci  •  *i.0e00E*8a9 

RMS  l/lo  •  ♦1.08BBE-BB2 

C 2  -  0. 0000E+000 

MEASURED  HEIGHT  ■  37.38 

FREQUENCY  •  30i 

CALCULATED  DEN9ITY  •  *3 . 9.26E-003 

•  8F  SAnPLES  •  Ctl 

MEAN  UOLTACE  •  «1.137CE*M( 

Figure  B.3  Typical  Fluctuations  in  Dust  Probe  Output 

Errors  in  the  measured  C  cannot  be  determined  by  repeated  measurements  made 
with  the  same  hopper  opening,  since  this  does  not  separate  different  realiza¬ 
tions  from  variations  in  the  measurement  of  C.  Instead,  the  uncertainties  of 
the  components  of  the  measurement  of  C  were  analyzed  as  follows.  The  uncer¬ 
tainty  in  C  is  given  by 


dC  3  Ki  “"l2  *  (§t  dV>2l‘/2 

1  /9 

dC  =  [(i  dm)2  +  (—  dV)2] 

V2 

r  ■  Ks2)2  +  (r-)2]172 


(B.12) 
(B . 13) 
(B.14) 
( B .  15 ) 
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dV/V  is  easy  to  calculate,  since  the  volume  of  the  tube  was  determined  by 
V=  ur2h,  where  r  =  7.26  cm  ±  2%,  h  =  58.7  cm  ±  1% 


■  2( . 02 ) 
■  7.26 


)2  + 


.01  f  1/2 

5S77j  J 


1% 


(B . 16 ) 


There  were  three  contributors  to  dm/m.  First,  there  was  about  1/8  inch  (3  mm) 
play  in  blades  relative  to  one  another,  thus  in  an  extreme  case  the  top  blade 
could  lead  or  lag  the  bottom  blade  by  1/8  inch  for  the  aproximate  100  ms  clo¬ 
sure  time,  causing  a  .36  in2  (2.3  cm2)  area  difference  to  be  exposed.  For  a 
density  C  flowing  at  a  velocity  V  (assumed  previously  to  be  30  cm/s)  this  A 
area  would  cause  a  A  mass  given  by  Am  =  Amt  =  C  VAAt,  or  since  C  =  mass/ 
vol ume. 


Am 


_  mVAAt 
"  ”701“  * 


(B.17) 


Therefore,  dm/m  can  be  approximated  by 


dm  _  Vel 
m  Vo T 


closure 


_  30  cm/s  x  2.3  cm2  x  0.1  s  =  7  x  10-4 


9727  cm3 


(B. 18) 


which  is  negligible.  The  second  contributor  to  dm/m  is  weighing  innaccura- 
cies.  The  weight  could  be  measured  to  ±  .02  grams.  Taking  the  worst  case,  a 
measurement  of  C  =  1  x  10"4  g/cc  with  the  tube  volume  of  9727  cm3  means  that 
.97  grams  are  measured  with  a  ±  2%  error.  The  third  contributor  to  dm/in  is 
spatial  average  non-uniformities  in  the  flow  during  the  sampling  time.  These 
were  not  directly  measure,  but  can  be  assumed  to  be  of  the  same  order  of  mag¬ 
nitude  as  the  temporal  variations  observed  in  the  voltage  signal,  i.e.,  ± 
10%.  Assuming  these  errors  in  dm/m  are  additive,  dm/m  =  12%.  Therefore 


=  ( ( . 12 ) 2  +  (.01)2 ) 1/2  =  12% 


(B.19) 


The  path  length,  h,  was  measured  with  dial  calipers  and  so  can  be  considered 
accurate  to  .01  inch  or  1%  for  the  short  path  length. 
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Uncertainties  in  Extinction  Curves  Calculated  From  Soil  Analysis 


The  uncertainties  in  the  extinction  term  are  given  by 

d32P 


dd32  2  ,  2  1/2 

t(T7T>  *  0 1 


(B.20) 


Since  values  for  d32  and  p  were  obtained  from  an  independent  soils  lab  and 
from  suppliers  data  sheets,  it  is  difficult  to  assess  the  uncertainties  asso¬ 
ciated  with  them.  One  way  to  estimate  the  magnitudes  is  to  look  at  the  dif¬ 
ferences  between  results  obtained  at  the  start  of  the  program  and  results 
obtained  at  the  end  of  the  program.  The  particle  specific  gravity  of  the  soil 
sample  was  given  as  2.57  before  the  test  and  2.59  after  the  test,  or  a  ~  1% 
change.  (No  pretest  data  on  the  Ottawa  sand  specific  gravity  was  available.) 
An  estimate  of  the  uncertainty  in  particle  density  is  then 


=  1Z 
P 


(B.21) 


The  value  of  d32  calculated  from  the  suppliers  data  sheet  for  Ottawa  F-140 
(Appendix  A)  was 


d32  =  .0078  cm 


(B .22 ) 


and  the  value  calculated  from  the  independent  laboratory  analysis  of  this 
Ottawa  sand  was  d32  =  0.0110  cm.  Since  the  Ottawa  sand  sent  to  the  laboratory 
was  from  a  virgin  (not  used  in  the  calibration)  sample  of  the  Ottawa  sand, 
this  will  be  one  estimate  of  the  uncertainty  in  d32 

dtd32)  _  .oil  -  .0078  cm  _ 

“7 - .UU94  cm - 3555  (B-23) 

d  3  2 

The  pre-test  and  post-test  results  obtained  for  the  Pre  Direct  Course  soil 
samples  were  d32  =  .0190  cm  and  d32  =  .0213  cm  respectively.  Then  a  second 
estimate  of  the  uncertainty  in  d32  is 


B-6 


The  uncertainty  in  the  term  6^2  can  &e  considered  to  lie  between  11%  and  35%, 


and  will  be  estimated  as  20%. 
. 0 1 2 ) 1/2  -  20%. 


The  error  in  the  term 


1 


d3  2  p 


is  the,  (.22  + 


APPENDIX  C.  MEASUREMENT  STATION  BLAST  LOADS 

The  blast  loading  on  the  measurement  station  was  modeled  on  the  information 
presented  by  Brode9.  Two  cases  were  investigated,  a  100  psi  overpressure 
ideal  shock  and  a  100  psi  overpressure  precursor  shock.  In  the  ideal  shock 
case,  the  pressure  jumps  to  the  peak  pressure  within  a  fraction  of  a  milli¬ 
second.  The  decay  of  the  pressure  pulse  is  initially  dominated  by  the  passage 
of  the  shock  front  itself,  then  the  further  decay  is  dictated  by  the  rate  of 
pressure  decrease  in  the  shock  wave  interior  (Figure  C.l). 


Time 


Figure  C.l  Fast-Rise  or  Ideal  Shock 

In  the  precursor  shock  case,  the  pressure  pulse  has  a  slower  rise  to  peak  and 
a  more  irregular  decay  than  that  exhibited  by  the  ideal  shock  (Figure  C.2). 


Time 


Figure  C.2  Slow  or  Finite  Rise-Time  Shock 


C-l 


I 


For  certain  diffraction  type  targets  such  slow  rising  pressures  can  reduce  the 
damage  potential.  For  drag  type  targets  the  damage  may  actually  increase 
because  higher  dynamic  pressures  develop.  The  pressure  histories  for  the  100 
psi  overpressure  ideal  shock  and  precursor  shock  were  calculated  to  have  the 
form  illustrated  below  (Figure  C.3) 


Tint  (ms)  TIME  (ms) 

Figure  C.3  Front-Face  Pressure  Versus  Time 

In  the  present  analysis,  the  probe  stings  and  rake  blades  were  analyzed  sepa¬ 
rately.  In  each  case  it  was  determined  whether  the  component  responded  to 
diffraction  loads,  drag  loads,  or  both,  and  then  the  stresses  and  deflections 
were  calculated. 
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The  dynamic  response  of  the  probe  stings  and  rake  blades  to  blast  loading  was 
approximated  by  an  undamped  single  degree  of  freedom  system  subjected  to  a 
transient  exciting  force.  Figure  C.4  represents  response  spectra  of  such 
undamped  oscillators  for  the  excitations  used  to  characterize  the  blast  loads. 
In  these  figures,  the  ratio  t^/x  is  the  ratio  of  the  characteristic  time  of 
the  exciting  force  divided  by  the  natural  period  of  the  oscillator.  The  term 
(xk/F0)max  is  a  nondimensional  number  which  is  a  measure  of  the  dynamic  effect 
over  that  of  a  statically  applied  load. 


r~  'i  — *4 


Figure  C.4  Response  Spectra  of  Representative  Excitations 


The  characteristic  time  of  the  reflected  pressure  shock  was  determined  by  the 
fact  that  the  reflected  pressure  shock  is  reduced  to  the  stagnation  pressure 
at  a  time  after  shock  arrival  approximately  equal  to  3S/U  where  S  equals  the 
height  or  half  width  of  the  structure  (whichever  is  less)  and  U  is  the  shock 
speed.  For  100  psi  overpressures  acting  on  the  measurement  station,  this  time 


is  on  the  order  of 


i*  ift. 


t1  =  1  x  10  "r  sec 


IC.1J 


Thus,  the  initial  pulse  of  the  ideal  shock  load  was  modeled  by  Figure  C.4(a) 

The  "tail"  of  the  ideal  shock  load  and  the  precursor  shock  have  the  same 
characteristic  time 


tj  =  4  x  10  sec 


(C.2) 


Ths  loading  was  modeled  by  Figure  C3 . 4 ( a )  if  x  ~  tj,  or  by  Figure  C.4(b)  if 
x  >>  tj  or  x  «  t}.  In  general,  if  the  component  was  found  to  respond  to  the 
blast  load  the  dynamic  effect  was  automatically  assumed  to  be  an  order  of  2 
greater  than  the  effect  due  to  a  statically  applied  load  in  order  to  be  on  the 


safe  side. 


PROBE  STINGS 


The  advantages  of  using  a  linearly  tapered  sting  as  opposed  to  a  uniform  dia¬ 
meter  sting  are  illustrated  as  follows.  The  tip  deflections  and  tip  angles  of 
a  tapered  cantilever  beam  relative  to  a  uniform  cantilever  beam  of  Ia  (Figure 
C.5)  are  given  for  the  case  of  an  end  load  by 


(C.3) 
(C.4) 

and  for  a  distributed  load  by 


y;  =  -354  *a 

9'  =  .402  8, 
a  d 


ya  =  .328  ya  (C.5) 

9,  =  .354  Ga  (C.6) 

d  a 


Therefore,  tapering  the  beam  provides  more  than  a  factor  of  2  decrease  in  the 
tip  deflection  and  angle.  In  addition,  this  removes  stress  raisors  and 
reduces  the  maximum  bending  stress  aTOX  =  Mc/I  at  the  base. 


Figure  C.5  Tapered  Sting 


Therefore,  the  sting  will  be  designed  with  a  linear  taper  from  an  outer 
diameter  of  3/4"  at  the  tip  to  an  outer  diameter  of  1"  at  the  base. 

The  need  for  a  symmetric  noise  piece  is  illustrated  as  follows.  If  the  probe 
has  a  knife  edge,  it  is  necessary  to  look  at  the  aerodynamic  tip  lift  due  to 
the  dynamic  pressure  occurring  in  ideal  shock.  This  lift  will  cause  a  normal 
force  at  the  tip  given  by 

Fn  =  4q  CD  A  =  500  lbs.  (C.7) 

Using  flexure  formula,  the  maximum  bending  moment  due  to  this  force  is 

Mb  =  FnA  =  500  lb  x  8  in  =  4,000  in  lb  (C.8) 

which  results  in  a  fiber  stress  in  the  sting  given  by 

a  --  =  14-0TO  1"  1b)  I-.5-.!"1  »  4.3  x  10-  lb/in2  (C.9) 

mx  ^  .046  in- 

This  is  unacceptable,  since  the  yield  strength  of  steel  is  33  x  103  lb/in2. 
This  points  out  the  need  for  a  symetric  nose  piece. 
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Elastic  Stability 


The  elastic  stability  of  the  stings  was  investigated  by  finding  the  critical 
buckling  load.  As  a  conservative  estimate,  a  uniform  beam  with  the  tip  outer 
diameter  was  used  rather  than  the  tapered  beam.  The  relevant  parameters  are 
then 


1/2"  1.0.  steel  tube  (C.10) 

3/4"  O.D.  (C.ll) 

I  =  (f)(Rj  "  Rj)  =  (J)(.3754  -  .254  )  =  1.24  x  10“2  in4  (C.12) 

A  =  it  (Rj  -  r2)  =  m  ( . 3752  -  . 2 52 )  =  2.45  x  10"1  in2  (C.13) 

E  =  30  x  106  lb/in2  (C . 14 ) 

=  0.25  (C.  15) 

which,  when  used  in  the  critical  load  formula 


=  5.8  x  104  lbs  (C. 18) 

Therefore,  each  sting  could  support  5.8  x  104  lbs  of  static  purely  axial 
loading  without  buckling. 
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Natural  Period  of  the  Stings 


The  transverse  and  axial  natural  periods  of  the  sting  were  found  by  modeling 
it  as  a  uniform  cantilever  beam  of  annular  cross  section  with  a  weight  at  the 
tip.  The  natural  period  of  the  tapered  beam  (Figure  C.5)  was  assumed  to  lie 
between  that  of  the  uniform  beam  with  the  tips  cross-sectional  area  A,  and 

a 

moment  of  inertia  Ia,  and  that  of  a  beam  with  the  base  crossectional  area  Ab 
and  moments  of  inertia  Ib. 

The  relevant  properties  for  the  sting  models  of  the  two  different  diameter 
beams  were  found  to  be 


0 

a)  r  =  .75" 

a 

0  b) 

rb  = 

E 

=  30  x  106  lb/in2 

E  = 

30  x  106  lb/in2 

(C.19) 

=  .0125  in1* 

h  = 

.035  in4 

(C.20) 

Aa 

=  .245  in2 

Ab  “ 

.59  in3 

(C.21) 

wa 

=  6.9  x  10-2  lb/in 

u>b  = 

.17  lb/in 

(C.22) 

g 

=  386  in/s2 

g  = 

386  in/s2 

(C.23) 

The  tip  weight  was  approximated  as  the 

weight  of 

half  of  a  steel  cube  with  .75 

in  sides,  giving 

W  =  9.5  x  10"2  lb. 

The  transverse  natural  frequency  of  a 

uniform  cantilever  with  a  tip 

weight  and 

distributed  weight  is  given  by 

f  =  . 

1.732  Elg 

(C  .24 ) 

1 

WA3  +  .236  u>A4 

which  gives  the 

transverse  natural  frequencies  of 

the  two  beams  as 

f  ~ 

A 

310  Hz  f. 

~  435  Hz 

(C.25) 

Since  the  response  to  fast  risetime  transient  excitation  is  the  important 
property,  the  estimate  will  be  weighted  toward  the  higher  frequency,  and 


considered  to  be 


f  ~  400  Hz  (transverse  excitation) 


(C.26) 


The  axial  natural  frequency  is  determined  by 


rl  "  TH 


WA  +  (dA2/ 3 


(C.27 ) 


which  gives  the  axial  natural  frequency  of  the  two  beam  models  as 


f  =  5,700  Hz 

a 


fb  =  6,300  Hz 


(C .28) 


so  the  axial  natural  frequency  will  be  considered  to  be 


fn  ~  6,000  Hz  (axial  excitation) 


(C.29) 


Thus,  the  natural  periods  for  the  axial  and  transverse  vibrations  of  the 
tapered  s*’ng  are  taken  to  be 


Taxial  =  *  10  4  s 


(C.30) 


^transverse 


=  2.5  x  10"3 


(C.31) 


Sting  Stresses,  Tip  Deflections  and  Tip  Angle 


The  ideal  shock  diffraction  load  (PR)  has  the  character!’ Stic  time 


=  1  x  10  s 


C-9 


and  since 


(C.33) 


(C . 34 ) 


the  axial  excitation  causes  a  dynamic  response  that  is  approximately  twice  as 
great  as  the  static  response.  The  reflected  pressure  will  cause  an  axial 
transient  loading  of 


F,  xial  3  pref  AS  cos  45  +  Pref  AS  cos  45  (C.35) 

Pretl  =  430  psi  (C . 36 ) 

As  (.75)  (— g-— )  =  .4  in2  (C. 37) 

Fa x i a  1  =  255  lb  <C-38> 

which  will  result  in  tension  and  compression  stress  waves  with  a  maximum 

°  ~  2astatic  =  2  ^axial  Across  sectional  2  2,077  lb/in2  (C.39) 


This  is  acceptable,  since  the  yield  strength  of  steel  is  33  x  103  lb/in2. 


The  ideal  shock  lift  and  drag  loads  are  character! zed  by 


Vxtrans 


4  x  10~2  s 
2.5  x  10'Z  s 


16 


(C . 40 ) 


^l^axial 


4  x  10~2  s 
1.7  x  10~4  s 


400 


(C . 41 ) 


Therefore,  both  the  axial  and  transverse  excitation  must  be  examined. 


The  transverse  tip  load  due  to  aerodynamic  lift  should  be  negligible  for  a 
symmetric  nose  piece  since  the  angle  of  attack  is  zero. 


The  axial  loading  due  to  aerodynamic  drag  will  be 
Faxial  =  =  4  q  CD  A 

where 

CD  =  1.5 

Ap  =  (.75  in)2  =  .56  in2 
q  =  120  psi 

Faxial  =  403  lb 

which  will  cause  a  maximum  stress  in  the  sting  of 
a  -  2aS£a|..jc  -  3,291  lb/in 


(C .42 ) 
(C.43) 
(C.44) 
(C.45) 
(C.46) 

(C .  47 : 


which  is  acceptably  less  than  the  33  x  10  3  lb/in2  yield  strength  of  steel. 


The  precursor  drag  and  lift  loading,  as  illustrated  below,  will  cause  flexure 
of  the  sting.  This  flexure  is  due  to  the  combination  of  both  axial  and  trans¬ 
verse  loading  on  the  sting.  The  transverse  loading,  using  the  cos2  law,  is 


found  to  be 

fn  =  Tn  =  4  q  A/*  CD  cos20  (C.48) 
q  =  120  psi  (C.49) 
Ml  -  .75  i n 2 / i n  (C.50) 
CD  =  1.5  (C.51) 
9  =  70°  (C.52) 
fn  =  =  54  lbs/in  (C.53) 

The  axial  loading,  again  using  the  cos2  law,  is  found  to  be 

F  =  4  q  Ap  CD  cos20  (C. 54) 
6  =  20°  (C.55) 
Ap  =  (.75  in)2  =  .56  in2  (C.56) 
F  =  360  lb  (C.57) 


This  combined  loading  causes  a  tip  deflection  ya  and  a  tip  angle  9a  in  a 
uniform  cantilever  given  by 


—  (1  +  KA  tan  K l  - 

k2f 


cosKJr 


(C . 58) 


(C.59) 


where  K  =  /F/EI 


(C.60) 


C — 1 2 


I 


Therefore, 


ya  =  7.7  x  10“2  in 

ea  =  1.3  x  10-2  rad  =  7.3  x  lCf1 


(C .61 ) 
(C.62) 


The  tip  deflection  and  angle  of  a  tapered  beam,  as  previously  discussed,  are 
reduced  relative  to  those  of  a  uniform  beam  of  the  smaller  diameter.  Thus 


ya  =  .328  ya  (uniform) 
0,  =  .354  0a  (uniform) 


(C .63) 
(C.64) 


and  for  a  slowly  applied  load 


ya  =  2.5  x  10-2  in 
ea  =  4.6  x  10"3  ra< 


rad  =  2.6  x  10_1 


(C.65) 

(C.66) 


For  a  suddenly  applied  load  these  will  be  increased  by  a  factor  of  2,  and  the 
tip  deflection  and  angle  will  be 


ya  =  5  x  10“2  in 

0,  =  5.3  x  10_1  degress 


(C.67) 

(C.68) 


Dynamics  Technology 
DT-3338-06 


Effects  of  Tip  Motion  Due  to  Blast  Loading 


The  effect  of  tip  bending  in  the  plane  of  the  light  path  is  illustrated  by  the 
following  figure.  It  is  evident  from  this  figure  that  the  worst  case  occurs 
when  either  both  mirrors  tilt  "in",  or  both  mirrors  tilt  "out"  relative  to  the 
plane  of  the  light  path  (in  other  words,  if  both  stings  bend  towards  each 
other  or  apart  from  each  other).  If  both  stings  bend  equally  in  this  manner, 
the  light  reflected  off  the  second  mirror  is  deflected  by  an  angle  Op  equal  to 
4  times  the  tip  angle.  When  both  stings  bend  in  the  same  direction  the 
angular  deflections  of  the  mirrors  tend  to  cancel  each  other. 


Mirror  at  Mirror  at 


Figure  C.6  Effect  of  Tip  Bending  in  the  Light  Path  Plane. 


<*  <*  syy. 


A  different  situation  holds  for  bending  out  of  the  plane  of  the  light  path. 
The  mirror  plane  is  at  45°  to  the  plane  which  bends,  so  if  the  probe  tip  bends 
up  at  an  angle  9^,  the  light  beam  will  deflect  9^  down  in  the  vertical  plane 
normal  to  the  sting.  If  the  other  probe  tip  also  moves  an  angle  0^,  the  end 
result  will  be  to  preserve  the  downward  deflection  of  the  beam  by  9^  as  it 
reflects  off  the  second  mirror,  as  shown. 

Using  the  tip  deflections  and  angles  calculated  in  the  previous  section,  the 
worst  case  effects  of  the  precursor  blast  loading  on  the  light  path  can  be 
found.  The  uniform  drag  load  from  below  causes  both  stings  to  bend  upward, 
deflecting  the  beam  downwards  by  .53  degrees  across  the  interrogation  path. 
This  will  cause  a  4.6  x  10-2  cm  deflection  at  the  second  mirror.  The  second 
mirror  preserves  the  0  =  0.53  declination  angle,  so  the  beam  will  be  deflected 
downwards  by  an  additional  0.19  cm  as  it  travels  down  the  receiving  sting. 
When  the  upward  deflection  of  the  tip  is  taken  into  account,  as  illustrated, 
this  results  in  a  .11  cm  deflection  at  the  lens. 


ya -deflect  ion 


Figure  C.7  Effect  of  Tip  Bending  Out  of 
the  Light  Path  Plane. 

This  deflection  can  easily  be  accomodated  by  an  aperture  that  is  1.5  mm 
smaller  in  diameter  than  the  beam  diameter  at  the  lens. 

In  addition  to  the  deflection  of  the  light,  the  effect  of  off  axis  light  pass¬ 
ing  through  the  lens  is  to  shift  the  focal  point  by  a  distance 


S  =  f  tane 


(C.69) 


where  f  is  the  focal  length  of  the  collector  lens. 


Thus,  with  a  20  mm  focal  length  lens,  the  maximum  calculated  tip  motion  would 
cause  the  focused  light  to  move  .019  cm  down  the  face  of  the  receiving  fiber. 
This  movement  can  be  accomodated  by  a  receiving  fiber  that  is  either  much 
larger  than  the  focused  spot  (so  that  this  movement  does  not  place  any  of  the 
spot  off  the  fiber  face)  or  much  smaller  than  the  spot  size  (so  that  this 
movement  does  not  cause  the  edge  of  the  spot  to  cross  the  fiber  face).  From 
system  power  considerations,  the  larger  fiber  is  more  desirable. 

In  summary,  the  deflections  of  the  sting  due  to  precursor  blast  loads  (the 
worst  case)  can  be  accomodated  without  appreciable  signal  degradation  if  the 
aperture  on  the  receiving  optics,  and  the  receiving  fiber  are  sized  as  pre¬ 
viously  described.  This  does  not  represent  a  significant  departure  from  the 
probe  configuration  used  in  the  laboratory  tests.  Thus,  based  on  this  simpli¬ 
fied  analysis,  the  sting  sizing  is  adequate  for  precursor  blast  loads. 
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RAKE  BLADES 


Loading 

The  transverse  natural  period  of  the  rake  blade  about  its  minor  axis  is  x  = 
6.0  x  10"-  s,  so  it  will  respond  weakly  to  the  initial  shock,  and  strongly  to 
the  dynamic  pressure  loading  in  the  direction  of  the  flow.  The  transverse 
natural  period  of  the  rake  blade  about  its  major  axis  was  calculated  as 

-c  =  3.5  x  10~3  sec  (C.70) 


Therefore,  its  response  to  shock,  as  determined  by 
1  -  1x10  s  =  2.9  x  10‘2 


x  3.5  x  10_J  s 


(C.71) 


is  negligible.  Its  response  to  the  dynamic  pressure  loading,  characterized  by 


2i  =  4  x  10~2  s 
T  3.5  x  10-3  s 


=  11 


(C.72) 


will  be  significant. 


First,  looking  at  the  flexure  and  stress  due 
the  flow,  the  maximum  load  on  the  frontal  area 
precursor  shock,  and  has  a  value 


to  loading  in  the  direction  of 
of  the  rake  blade  occurs  in  the 


F  =  2Pr  As  sin(a)  +  Ap  Ap  +  CDAp  q 


(C.73) 


where  Pr  =  450  p si 

Ap  =  26  in2 
q  =  120  psi 
Ap  =  100  psi 
A$  =  26  in 
a  =  30° 

CD  =  2 


(C.74) 
(C.75) 
(C . 76) 
(C.77) 
(C.78) 
(C.79) 
(C. 


F  =  2.2  x  104  lbs 


(C . 81 ) 


; 


•*> 

i 

i 

* 


r 


r  '■ 


The  distributed  load  on  the  front  of  the  blade  is 

f  =  =  1.7  x  103  Ib/in  (C.82) 

If  a  worst  case  is  assumed,  in  which  this  distributed  load  continues  to  the 
blade  base,  the  top  of  the  rake  blade  will  deflect  backwards  by 

ya  =  ITT  (C.83) 

y  =  ILL x_  1?3  lb/ln)  (1,5  m)4 - =  l  s  x  10-3  in>  (C .84 ) 

(8)  (30  x  106  lb/in2)  (236  in4) 

If  this  is  increased  by  a  factor  of  2  because  of  sudden  loading,  the 
deflection  will  be 

ya  =  3.0  x  10~3  in.  (C.85) 

which  is  considered  negligible. 

The  maximum  fiber  stress  occurs  at  the  front  edge  of  the  blade,  and  is  given 
by 


a 


_  Me  _  f!2  c 

"  T~  "  ~T 


(C.86) 


a  =  U-7  x  103  Win)  (15_in)2  ( 6AJv±_  =  5>2  x  1q3  lb/in2  (c>87) 

(236  in4) 

If  this  is  increased  by  a  factor  of  2  due  to  the  dynamic  loading,  the  maximum 
tensile  stress  will  be  1.1  x  104  lb/in2,  which  is  sufficiently  less  than 
steels  yield  strength  in  tension  of  3.3  x  104  lb/in2. 


C — 1 3 


z  =  (2.2  x  104  Ib/in)  (15  in)4 

(8)  (30  x  106  Ib/in2)  (7.1  in4) 


=  0.65  inches 


(C  .95 ) 


which  is  unacceptable.  The  situation  can  be  rectified  by  bolting  cross  mem¬ 
bers  between  the  tops  of  the  blades.  These  cross  members  will  need  to  support 
half  the  total  load  on  the  blade  sides,  as  illustrated  below. 


Figure  C.9  Connected  Blade  Rakes 

or  F  =  1.1  x  104  lbs.  Since  the  yield  strength  of  steel  intension  is  3.3  x 
104  lb/in2,  the  cross  members  must  have  a  cross  sectional  rea  of  at  least 
.33  in2.  A  safety  factor  of  three  puts  the  requirement  at  1  in2,  which  can 
easily  be  accomplished  with  two  1/2  inch  by  1  inch  bars.  The  shear  on  the 
cross  member  bolts  will  be  5.5  x  103  lbs  per  cross  member,  which  can  be 
accomodated  by  two  5/8"  diameter  bolts. 

As  a  conservative  estimate,  if  the  precursor  flow  strikes  rake  blade  at  10° 
off  axis,  a  force  normal  to  the  side  of  the  blade  will  develop,  given  by 


F/R  =  4  q  CL  A 


(C .96) 


%ax  =  120  Psi 
Cj_  ~  0.6 

A  =  156  in2 


(C .97) 
(C .98) 


(C .99) 


E 


F  =  4  (120  psi)  (0.6)  (156  in2)  =  4.5  x  105  lbs 


(C.100) 


I  V 
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If  a  worst  case  is  assumed  in  which  this  distributed  load  continues  to  the 
blade  base,  the  static  tip  deflection  is  given  by 


y  =ZUl 

ya  8  El 


(C.102) 


=  (3.5  x  103  lb/in)  (15  in)4  =  1>Q  x  1Q-1  in 
(8) ( 30  x  106  lb/in2)(7.1  in4) 


(C. 103) 


With  a  dynamic  factor  of  2;  the  top  of  the  blade  will  move  ya  =  2.0  x  10_1 
inches. 


This  deflection  is  reduced  by  the  cross  members  at  the  top  of  the  blades.  The 
effect  of  this  modification  can  be  approximated  as  follows.  The  second  blade 
is  "shaded"  by  the  first  blade,  as  illustrated. 


Figure  C.10  Shading  of  One  Blade  by  the  Other 
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so  it  will  be  assumed  that  no  aerodynamic  forces  act  on  it.  When  the  top  of 
the  first  blade  moves  ya,  the  second  blade  resists  this  with  a  force 


F  =  K  y  where  K  = 

a  A3 


K  =  ( 3 )  (30  x  106_,1b/in2)  (7._l_iV  )_  _  (L9  x  105  i b/in ) 

(15  in}3 

The  situation  is  then  modeled  (as  shown  in  Figure  C.ll)  by 


(C . 104) 


(C.105) 


Figure  C.ll  Blade  Interaction 


using  the  principle  of  superposition,  the  deflection  due  to  the  distributed 


-f  A4 
ya  =  8EI 


and  the  deflection  due  to  the  restoring  action  of  the  second  blade 


"  ■ F  * 


(C.106) 


(C. 107 ) 


combine  to  give  a  total  deflection  of 


v 


(C.109) 


Solving  for  ya 


fa 

"W 


-]/[!  + 


K  A3 

EFT 


■1 


y 


a 


1.0  x  10  _ 

“  1  +  T7004  " 


.05  inches 


If  the  dynamic  factor  of  2  is  included 


(C.110) 


y  =  0.10  inches  (C . Ill ) 

a 

Thus,  the  cross  member  reduces  the  deflection  bu  a  factor  of  2.  In  reality,  a 
cross  member  would  have  even  more  effect  since  it  would  act  to  form  a  frame, 
and  thus  also  introduce  restoring  moments  at  the  blade  tips. 
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BOLT  SIZING 


Bolt  Shear  for  Precursor  Shock  with  Dust  Loading 

The  shear  on  the  bolts  holding  the  blades  to  the  base  plate  is  found  as 
follows.  The  total  axial  force  on  the  rake  was  calculated  in  the  discussion 
of  rake  blade  loading  as 

F  =  2.5  x  104  lbs  (C. 112 ) 

and  the  lift  force  on  the  blade  sides  was  found  to  be 

F  =  2.2  x  104  lbs  (C.113) 

and  thus  the  shear  at  blade  base  is 

V  =  3.3  x  104  lbs  (C.114) 

The  present  design  uses  (3)  1-1/4"  bolts  and  (2)  3/8"  bolts  to  secure  the  rake 
blade  to  the  base  plate.  The  maximum  allowable  shear  per  bolt  is  given  by  the 
bolt  area  x  10  Kips  so  the  maximum  allowable  shear  with  this  configuration  is 

V  =  (3)  (9.94  x  103  lbs)  +  (2)  (1.1  x  103  lbs)  =  3.2  x  104  lbs  (C.115) 

This  provides  33%  factor  of  safety  for  a  slowly  applied  load,  but  does  not 

include  the  factor  of  2  for  sudden  loading.  This  factor  of  2  can  be  obtained 
by  adding  a  "backbar"  with  (2)  1-1/4"  bolts  per  blade  to  bring  the  total 
number  of  1-1/4"  bolts  per  blade  to  (5).  Thus,  the  maximum  allowable  shear  is 
increased  to 

V  =  (5)  (9.94  x  103  lbs)  +  (2)  (1.1  x  103  lbs)  =  5.2  x  104  lbs  (C.116) 

which  provides  the  factor  of  2  necessary  for  sudden  loading. 
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BOLT  TENSILE  STRESS 


The  tensile  stress  on  the  bolts  that  hold  the  base  plate  to  the  floor  (or  the 
concrete  pad_  can  be  found  as  follows.  If  the  axial  force  on  the  blades  found 
in  the  previous  section  is  considered  a  distributed  load, 

f  =  =  1.8  x  103  Ib/in  (C.117) 

and  as  a  worst  case  this  distributed  load  acts  over  the  entire  blade,  the 
moment  is 


M  =  =  (1.8  x  103  Ib/in)  (15  in )2/2  =  2.1  x  105  lb/in  (C.118) 


The  maximum  tensile  fiber  stress  in  the  rake  blade  occurs  at  the  leading  edge- 
of  the  rake  base,  and  has  a  magnitude 


„  _  Me 

3max  '  T~ 


(C.119) 


o  =  (2.1  x  105  lb  ,in,M6.4  in)  •  6>Q  x  1q3  lb/in2 
(236  in4) 


(C. 120) 


which  is  approximately  a  factor  of  5  below  steel's  yield  strength  of  3.3  x  104 
lb/in2  in  tension. 

The  tensile  stress  in  the  base  plate  bolts  can  be  approximated  by  the  taking 
the  moment  due  to  the  actual  force  on  the  rake  times  the  lever  arm  of  the  base 
At  a  point  7  in.  forward  of  the  rake  neutral  axis  the  upward  force  is 

F=( 6.0  x  103  lb-in)/ (7  in)  =  8.6  x  102  lbs  (C.121) 


» 
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Thus,  if  only  (1)  1-1/4"  steel  bolt  (.994  in2)  per  blade  is  placed  at  this 
point,  the  axial  stress  is 

a  =  liAjLi.9.1 1  bs  =  8>7  x  i03  lb/in2  (C.122) 

0.994  in2 


which  is  well  below  the  yield  stress  of  33  x  103  lbs/in2  for  steel. 


STRESSES  IN  THE  OPTICAL  WINDOW 


The  stresses  in  the  windows  are  approximated  by  the  stresses  in  a  simply 
supported  uniform  circular  plate  with  a  uniformly  distributed  load,  i.e.. 


Max  a 


(C . 123) 


where  t  is  the  thickness  (~  .125  inches) 
Me  is  moment  at  the  center  given  by 


M  =  (3+v) 


(C.  124) 


where 


q  is  the  uniformly  distributed  load  (100  psi ) 
a  is  the  diameter  (~  .6  in.) 
v  is  Poisson's  ratio  (~  0.28) 


M  =  (100  Psi)(.36  in2)_  ^  ^  =  ^  jp-lb 


(C.125) 


and  the  maximum  bending  stress  is 


Max  a  =  (6)  7-4  in~1b/in  s  2,834  psi 
(.125  in)2 


(C.126) 


Since  glass  has  a  compressive  working  strength  >  100,000  psi  and  tensile 
strength  >  25,000  psi,  this  is  considered  acceptable. 
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Dep  Ch  of  Staff  for  Ops  &  Plans 
ATTN:  DAMO-NCZ 

Dep  Ch  of  Staff  for  Rsch  Dev  &  Acq 
ATTN:  DAMA-CSS-N 

Harry  Diamond  Laboratories 

ATTN:  DELHD-NW-RH,  22800 

US  Army  Ballistic  Research  Labs 

ATTN:  DRDAR-BL,  R.  Eichelberger 

US  Army  Material  &  Median ics  Rsch  Ctr 
ATTN:  DRXMR-HH 


DEPARTMENT  OF  THE  ARMY  ^Continued) 

US  Army  Material  Command 
ATTN:  DRCDE-D 

US  Army  White  Sands  Missile  Range 

ATTN:  STEWS-NR-CF,  R.  Valencia 
ATTN:  STEWS-TE-LD,  J.  Blaine 

USA  Missile  Command 

ATTN:  DRSMI-RH 

DEPARTMENT  OF  THE  NAVY 

Naval  Ordnance  Station 

ATTN:  J.  Miller,  Code  5422  I 

Naval  Research  Laboratory 

ATTN:  Code  2627,  Tech  Lib 

Naval  Sea  Systems  Command 
ATTN:  SEA-0351 

Naval  Surface  Weapons  Center 
ATTN:  Code  K06 
2  cy  ATTN:  Code  K82 

Ofc  of  the  Deputy  Chief  of  Naval  Ops 

ATTN:  NOP  654,  Strat  Eval  &  Anal  Br 

Strategic  Systems  Project  Office 
ATTN:  NSP-272 

DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Geophysics  Laboratory 
ATTN:  D.  McLeod 

Air  Force  Rocket  Propulsion  Lab 
ATTN:  LKCP 

Air  Force  Wright  Aeronautical  Lab 


Air  Force  Systems  Command 
ATTN:  OLW 

Air  Force  Weapons  Laboratory 
ATTN:  SUL 

Air  Force  Wright  Aeronautical  Lab 
ATTN:  FIBC 
ATTN:  FI MG 

Air  University  Library 
ATTN:  AUL-LSE 

Ballistic  Missile  Office 
ATTN:  ENMR 
ATTN:  ENSN 

Deputy  Chief  of  Staff,  Rsch,  Dev  &  Acq 
ATTN:  AF/RDQI 
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DEPARTMENT  OF  JHE_  AIR  FORCE  (Continued^ 

Deputy  Chief  of  Staff,  Rsch,  Dev  &  Acq 
ATTN:  AFRDS,  Space  Sys  &  C3  Dir 

Foreign  Technology  Division 
ATTN:  SDBG 

Strategic  Air  Command 

ATTN:  NRI/ST INFO 

Strategic  Air  Command 
ATTN:  XPFS 

Strategic  Air  Command 
ATTN:  XPQ 

161  ARG  !  izona  Ang 
2  cy  ATTN:  Lt  Col  Sherer 

OTHER  GOVERNMENT  AGENCIES 

Central  Intelligence  Agency 
ATTN:  OSWR/NED 


NATO 

NATO  School,  SHAPE 

ATTN:  US  Documents  Officer 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

University  of  California 

Lawrence  Livermore  National  Lab 
ATTN:  L-10,  H.  Kruger 
ATTN:  L-8,  P.  Chrzanowski 

Los  Alamos  National  Laboratory 
ATTN:  R.  Selden 

Sandia  National  Laboratories 
ATTN:  0.  Ri gal i 
ATTN:  Org  7112,  A.  Chabai 
ATTN:  R.  Clem 
ATTN:  0334,  J.  Struve 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 


Acurex  Corp 

•  >, 

ATTN: 

C.  Nardo 

ATTN: 

C.  Powars 

»/.“ 

ATTN: 

J.  Crenshaw 

a* 

ATTN: 

J.  Saperste 

u 

Aerojet  General  Corp 

y\ r: 

ATTN: 

R,  Jenkins 

\y- 

Aerospace  Corp 

•V  •_ 

ATTN: 

D.  Platus 

-V- 

ATTN: 

W.  Portcnier 

APTEK,  Inc 

ATTN: 

T .  Meagher 

AVCO  Systems 

Division 

ATTN: 

A.  Pal  lone 

1  ' 

ATTN: 

J.  Stevens 

ATTN: 

V.  Dicristina 

' : ' 

ATTN: 

W.  B coding 

I'*.*. 
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DEPARTMENT  OF  pEFENSEJlONTRACTORS  JContinuedJ 
Dynamics  Technology,  Inc 


2  cy 

ATTN: 

M.  Petach 

2  cy 

ATTN: 

C.  Dube 

2  cy 

ATTN: 

D.  Hove 

Fiber 

Materials,  Inc 

ATTN: 

M.  Subi 1 ia  , 

General  Electric  Co 

ATTN: 

A.  Garber 

ATTN: 

B.  Maguire 

ATTN: 

D.  Coceano 

General  Research  Corp 

ATTN: 

R.  Parisse 

ATTN: 

R.  Williams 

Institute  for 

Defense  Anal 

ATTN: 

Classified  L 

Kaman 

Sciences  Corp 

ATTN: 

F.  Shelton 

Kaman 

Sciences  Corp 

ATTN: 

E.  Conrad 

Kaman 

Tempo 

ATTN: 

DAS  I  AC 

Kaman 

Tempo 

ATTN: 

DAS  I  AC 

Lockheed  Missiles  &  Space 

ATTN: 

D.  Price 

ATTN: 

G.  Chrusciel 

ATTN: 

0.  Lee 

ATTN: 

P.  Schneider 

Martin 

Marietta  Corp 

ATTN: 

L.  Kinnaird 

McDonnell  Douglas  Corp 

ATTN: 

L.  Cohen,  MS 

2  cy 

ATTN: 

H.  Morris 

2  cy 

ATTN: 

J.  Mardacci 

McDonnell  Douglas  Corp 
ATTN:  M.  Potter 

National  Academy  of  Sciences 

ATTN:  Natl  Materials  Advisory  Board 

Pacific-Sierra  Research  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 

Pan  Am  World  Service,  Inc  2 

ATTN:  AiDC/Library  Doc,  TRF 

PDA  Engineering 

ATTN:  j.  Dunn 
3  cy  ATTN:  C.  Thacker 

Physics  International  Co 
ATTN:  J.  Shea 

Rand  Corp 

ATTN:  P.  Davis 
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ATTN: 

F.  Field 

ATTN: 

D.  Landaur 

ATTN: 

ATTN: 

P.  Haas 

P.  Rausch 

Spectron  Development  Labs, 

ATTN: 

R.  Ross 

ATTN: 

T.  Lee 

D 

ATTN: 

G.  Gurtman 

SRI  International 

ATTN:  D.  Curran 

DEPARTMENT  OF  DEFFNSE  CONTRACTORS  (Continued) 

Rand  Corp 

ATTN:  B.  Bennett 

R&D  Associates 


S-CUBEO 


Science  Applications  Inti  Corp 
ATTN:  J.  Warner 

Science  Applications  Inti  Corp 
ATTN:  R.  Potts 

Science  Applications  Inti  Corp 
ATTN:  J.  Cockayne 
ATTN:  W.  Layson 

Science  Applications  Inti  Corp 
ATTN:  A.  Martellucci 

Southern  Research  Institute 
ATTN:  C.  Pears 

Cal  span  Corp 

ATTN:  M.  Holden 


D E PARTMENT  OF  DEFENSE  CONTRACTORS  (Conti nued) 

Space  Data  Corp 

ATTN:  B.  Bollermann 

Sparta  Systems,  Inc 


ATTN:  G.  Abrahamson 

TRW  Electronics  &  Defense  Sector 
ATTN:  D.  Baer 
ATTN:  R.  Myer 
ATTN:  W.  Wood 

TRW  Electronics  &  Defense  Sector 
ATTN:  E.  Allen 
ATTN:  E.  Wong 
ATTN:  L.  Berger 
ATTN:  V.  Blankinship 
ATTN:  W.  Grabowski 
ATTN:  W.  Polich 

Weidl inger  Assoc.,  Consulting  Engrg 
ATTN:  T.  Deevy 
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